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MATHEMATICAL BIOLOGY OF DIVISION OF LABOR 
BETWEEN TWO INDIVIDUALS OR TWO 
SOCIAL GROUPS 


N. RASHEVSKY 


COMMITTEE ON MATHEMATICAL BIOLOGY 
THE UNIVERSITY OF CHICAGO 


The concept of the satisfaction function is applied to the situation 
in which two individuals may each produce two needed objects of satis- 
faction, or may each produce only one of the objects and then make a 
partial exchange. It is shown that with a logarithmic satisfaction func- 
tion there is no advantage in a division of labor, unless such a division 
materially increases the purely physical efficiency of production. This re- 
sult appears to be connected with the particular choice of the form of sat- 
isfaction function (logarithmic). While the problem has not been solved 
for other forms, it is made plausible that satisfaction functions which have 
an asymptote will lead to a different result. 

Next the case is studied in which division of labor occurs between two 
groups of individuals. It is shown that in this case the relative sizes of the 
two groups are determined from considerations of maximum satisfaction. 
Possible applications to problems of urbanization are suggested. 


In our two books (Rashevsky, 1948, referred to as HR; 1951, referred 
to as SB) we have used the concept of the satisfaction function to study 
problems of motivational interactions between two or more individuals. 
Other studies in that direction have been made by Anatol Rapoport 
(1947a, b, c). All the studies were limited to the case in which all individu- 
als engage in the same activity and produce, as a result of it, the same 
object of satisfaction. The parameters of the satisfaction function were 
considered as different for different individuals. 

In this paper we shall study a few simple instances of a different case— 
that in which two or more objects of satisfaction are involved. We shall 
restrict ourselves here to the case of two objects of satisfaction only. As 
before, we shall assume that the satisfaction function consists of two parts. 
One varies logarithmically with the amount of the object of satisfaction, 
and is positive for sufficiently large amounts; the other varies linearly with 
the amount of work necessary to produce the object of satisfaction, and 
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is always negative. That work is most appropriately measured by the ac- 
tual energy expenditure. However, we shall measure it here in terms of the 
time spent on work, which is equivalent to assuming that for different 
tasks the energy expenditures per unit time are the same. Though this is 
actually not so, this limitation in no way affects our results. 


i 


Let x be the amount of one object of satisfaction, and y of another. Let 
the amounts of time needed to produce x and y be v and w respectively, 
and let the amounts of objects of satisfaction be proportional to the times 
(HR, chap. xix; SB, chap. xv). Thus, denoting by a and a2 two constants, 


we have 
Y= OW 2 y = dw . (1) 


Let an individual be engaged in the production of the first object only. 
Then (HR, chap. xix; SB, chap. xv), denoting by A, and B, two other con- 
stants, we have for his satisfaction 


S’ = A, log a,9—B,0. (2) 
This has a maximum for 
een (3) 


Since the individual adjusts his work so as to maximize his satisfaction, 


the latter will be 
S’ =A, log Ba Ay. (4) 


If the same individual engages in the production of the second object of 
satisfaction, then, denoting by A: and B, two constants, we have 


S” = A» log Qow — Baw , (5) 
which has a maximum for 
— A 2 
and the maximum value S’’ is : 
ae A»yd2 
Sy, = As log —— Aa. . (7) 


The total time which an individual can devote to any work (or, for that 


matter, the total energy) is limited. Let us denote that maximum value 
Dygel oat 


Ai as tiugt < : 
Bt RB, <0: | (8) 
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then the individual can engage in the production of both objects of satis- 
faction, and maximize his satisfactions S’ and S"’ for both objects inde- 
pendently. He will have an amount A,a;/B, of the first and an amount 
A2d2/By of the second. 

If, however, 


B44 (9) 


Bs 


then the individual cannot maximize S’ and S” independently, because 
v and w are not independent. We now have 


v+w=T. (10) 


Eliminating w from expression (5) by means of (10), we now find for the 
total satisfaction 


es eS! =A; log 0, — Bi 7-+ As log G2 (T— 2), —Ba(T —i0).. (11) 


Now the individual can adjust his only variable v so as to maximize S by 
requiring 
——=0. (12) 


Introducing (11) into (12), we find 
A,(T — v) —Biyv (T— v) —Agvt+tB.y(T—v) =0, (13) 


which is quadratic in v, and leads to somewhat unwieldy expressions. For 
simplicity, we shall restrict ourselves to the case Bj = B, = B. In other 
words, we assume that the unpleasantness of both types of work is the 


same. We then find 
Ay 


eee 7 
A;+ A. 


Introducing this into (11), we find for the maximum value S* of S 


v= (14) 
a 
S*=A; log Me ala eaay AT = pip (15) 


If inequality (8) holds, then the maximum total satisfaction S,, is again 
simply given by the sum of expressions (4) and (7). Thus 


Sm = Ay log “34 As log As (Ay + An). (16) 


We now shall prove that - am | 
ig | SS Sy (17) 
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Introducing (15) and (16) into (17) and rearranging, we see that the 


latter is equivalent to 


A A 
(A+ As) log B+ (Art-As) < (Ar+As) log “2 FE“? +7. (18) 


After further rearrangements, inequality (18) may be brought into the 
form 


log <a! (19) 
Put 
1— eat =e. (20) 
But S* has a meaning only if (9) holds. This implies that 
; Oe <1: (21) 
Inequality (19) now becomes 
log. == ey <a aren (225) 


That this inequality always holds is seen from the fact that the expan- 
sion of log (1 — z) consists of only negative terms, the first being —z. 
This proves inequality (17), which is also intuitively evident. 

Now consider two individuals, each in the position described above, 
that is, each producing both objects of satisfaction under the condition 
expressed by inequality (9). However, while (9) holds, let 


= aie and = 5 Be (23) 

We may now ask whether the two individuals will gain anything by 
dividing their labor, that is, one of them producing only the first object 
of satisfaction, the other the second, and then sharing their results in an 
appropriate manner. 

Let the first individual produce only the first object of satisfaction in 
an amount %o, while spending the time vo; the second individual produce 
only the second object of satisfaction in the amount yo, spending the time 
wo. Then let them divide the results so that the first will have x, and 
of both objects, the second a and y,. Thus 


Xy + X= Xp 5 (24) 
Yitve= Vo. (25) 


Since the first individual works only to produce the first object of satisfac- 
tion in the amount 2, therefore, the negative part of his satisfaction func- 
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tion will always be Byy = Bxo/a;, regardless of the method of sharing. 
Similarly, the negative part of the satisfaction function of the second in- 
dividual will be Bwo = Byo/a:. The total satisfaction functions will be, 
respectively, 


S, = A, log %1+ Ao log y= Xo; (26) 
1 
B 

pes 4 log tat dalog y2——— e- (OP) 


How will the quantities «, a2, 1, and ye be determined? This problem of 
“trading” one object for another, between two individuals, when each 
tries to obtain the best advantage from the trading, has been solved for 
two individuals and two objects by G. Evans (1930, pp. 125-28). A more 
general case, for any number of individuals and any number of objects, 
has been treated by the present author (1952a). The reader is referred for 
details to the two sources mentioned above, where the proof of the follow- 
ing is given. 
We form the expressions 


A1==— } Yi=——; X2=3—; Y,=——; (28) 
1 


We then have in the present case 


X1(%1, 91) = V1 (%1, y1) (29) 
X2 (Xp — X1, Yo— M1) Vo (% — %1, Yo— M1)’ 
and 
Nise — Xi (41, Yo) (30) 
x1 — Xo V1 (4%, 1) © 


Equations (29) and (30), together with (24) and (25), determine the four 
unknown quantities, #1, %2, V1, and 4». 
From (26) and (27) we have 


As xX nae V Ao. V Ay (31) 


Introducing this into (29) we obtain, after rearrangements, 


eee (32) 
Hy, M1 


Introducing (31) into (30), we find 
ee yo ape aioe (33) 
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Equations (24), (25), (32), and (33) give 


WA eee | 
X71 A Ao 0%. | +A, 0° 
fais 1 (34) 
ae a 
Vi eA; ee 
Introducing these values into (26) and (27), we find 
= A1%9 AV ome oa ae 25 
SiS Ar logan gy aes ey - Xo 5 ( ) 
ay AoXq Ayo Bez 36 
Paper parr se Oko sare oe Ce 


The amounts a and yy are still left indefinite. Expressions (35) and (36) 
specify S; and S; for any given values of «» and yo. These will be deter- 
mined by the two individuals. If they act egoistically, or individualistical- 
ly, then the first individual will choose his x so as to maximize S;, while 
the second will choose such a yo as to maximize S2. In other words, %» and 
yo are in this case determined by 


Equations (37), together with (35) and (36), give 
Ajay r Ad 


Opa eens ae 


B ? BA Serena) (38) 


which are the same as would be obtained from (3) and (6), for the case of 


a single individual working under conditions (8). Introducing (38) into 
(35) and (36), we obtain 


Ajay A,Aoa 

Sim A Om ee ee ee eee . 

; 1 °°8 (Ay Ag) 8 Bae 
A,Aoa Aja 

Som = 4 loge ee poe ME Ne eas 

: 198 B (A; ADIL ee 


In order to see whether the two individuals gain anything from the di- 


vision of labor, we must compare Sim and S2m with S* as given by (15). 
We have, from (15) and (39), 


Sim —S* = (Ait A») log 214 (BE A). (41) 
Put 
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Because of (9) and (23) 
Oe lee 


Now (41) may be written 
Sim—S* = (A,+A2) log (1 — 2) +BTz. (43) 
Expanding log (1 — z), we find 


Sim — S* = [BT — (A14+ A2) |] 2-+ negative terms in 2. (44) 
Hence 
Sim—-S* <0. (45) 


Similarly, we prove that 
aoa ee Ue (46) 


Thus in the case of egoistic behavior, the division of labor does not lead 
to any advantage, with the assumptions made above. 

If both individuals exhibit altruistic or collectivistic behavior, they will 
determine x and yo from 


kort one on) OC Sle Sz) ay 
Introducing (35) and (36) into (47), we now find 
_ 2 Aid | weArh 
DA poeta pt dl aaa la (48) 


Introducing (48) into (35) and (36), we obtain corresponding expres- 
sions for Sim and Som, which now differ from (39) and (40) by (41 + Ae) 
log 2. Therefore, instead of (44), we now find 


Sin» = (Ai Aa) lop 2(1— 2) BT 2. (49) 
ieee, that is, if 
Ll 5 (50) 
Bees 


then Sim — S* > 0. It will, however, be positive even for values of z which 

are not too much smaller than 3, for then the positive term BYs in (49) 

will still prevail. Thus for altruistic behavior Si, > S* if (50) holds, and, 
sunilarly, Son, > S* if ee, 

2 

RB = 2 ie (Sal) 

We notice that at least one of the two inequalities (50) and (51) must 

hold, if (9) holds, with B; = B, = B. Hence, for altruistic behavior, divi- 

sion of labor leads to an advantage for at least one of the two individuals. 

The important question is whether the above conclusions are due to the 

specific choice of the logarithmic function for the satisfaction. That choice 


leads to the following peculiarity. 
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If an individual produces an object of satisfaction, but has to give part 
of it away, retaining only a fraction «1, then his satisfaction function is 
given by 


S = A log ax — x. (52) 
This has a maximum for 
Aa 
Pe helias Sie: 
x B ? ( ) 


which is independent of a. With such a satisfaction function, an individual 
who has to share the results of his activity does not increase the activity. 
This results in a certain loss of advantage in sharing. If we put 


S=A log (1+ ax) 7 (54) 
then the optimum value of x is 
t=, (SS) 


In this case an individual will work the less, the more he has to share. 
A similar situation is obtained if we put 


S=A Vato a. (56) 


Such a situation is in no way implausible. With other satisfaction func- 
tions, the situation may be different. For example, put 


Sa iments g. (37) 


The value of « which maximizes (57) is 


1 Ada 
wheel log ae (58) 
For a = 1, when the individual retains all he makes, the maximum 
value x’ of x is 
il Aa 


‘— — 
a log 3° 


A (59) 


In order that a maximum exists at all in both cases, it is necessary that 


B 
Aa>B and peer ira (60) 
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From (58) and (59) we see that « > «’ implies 
ra 
B 

For a = 0, the value of —log a is positive infinite; for a = 1, it is zero, 
and the slope of the —log a curve at a = 1 is equal to 1 (Fig. 1). 

The straight line (1 — a) log Aa/B is tangent to the line —log a if log 
(Aa/B) = 1, or Aa/B =e (Fig. 1, broken line). In that case inequality 
(61) does not hold, and hence x’> «x. If, however, log Aa/B > 1, or 


Aa 
pes (62) 


then in the interval between a = a, anda = 1 (Fig. 1, full line) inequality 
(61) holds, and x > x’. The stronger inequality (62), the smaller the value 
of a;. Hence, if (62) holds, an individual who retains only a fraction a of 


(1—a) log >—loga. (61) 


N 
SS 


‘ 


fa) ox, | —> x 


FIGuRE 1 


what he produces will produce somewhat more than if he retains every- 
thing, provided a does not fall below a critical value a). That critical value 
decreases with increasing “efficiency” a, and with decreasing difficulty B 
of the work. If a < a, then the individual produces less than if he retains 
everything. If a < B/Aa [cf. (60)], then the individual will not work at all. 

It would be interesting to investigate the problem of division of labor 
with a satisfaction function of the form (57). Unfortunately, this leads to 
quite unwieldly expressions and transcendental equations. We shall, 
therefore, not investigate this problem now. Instead we shall make the im- 
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portant remark that very plausible conditions which affect the satisfaction 
function indirectly may result in an advantage of a division of labor for 
both individuals, even for egoistic behavior. An attempt by an individual 
to engage in two different tasks is likely to result in a decrease of efficiency 
of each task. This may be due to the necessity of special training in the 
particular task, or to the waste of time involved in switching over from 
one task to another. As a result of this, the values a; and a2 in expressions 
(39) and (50), which correspond to a division of labor, are likely to be 
much larger than the values of those quantities in expression (15) for S*. 
We shall denote those latter values by aj and ay. 
Instead of (41) we now have 


Ay pan ay 
Sim —S* = (A,+A,) log art (BT — A,) + A1 log a* 
ep rites) 
+ A, log Te: 


2 


As we have seen above, the sum of the first two terms is always negative. 


If, however, 
a,>a* ; a,>as, (64) 


then, if the last two inequalities are sufficiently strong, the last two posi- 
tive terms in (63) will prevail over the negative terms, and Sim — S* will 
be positive. The same holds for S:,, — S*. A division of labor may thus be 
of advantage to both individuals for purely physical, rather than psycho- 
biological, reasons. 


Ty; 


We now shall discuss the problem of division of labor not between two 
individuals, but between two social groups. 


Let the total population V9 of a society be divided into two groups, con- 
taining V; and NV, individuals, correspondingly. Let the individuals of the 
first group engage in the production of the first object of satisfaction, pro- 
ducing a total amount 2 for the group. Similarly, let all individuals of the 
second group engage in the production of the second object of satisfaction 
only, producing a total of yo. After a trading or exchange, such as dis- 
cussed above, the groups will possess amounts 21, 22, yi, and yo. For sim- 
plicity, we shall assume that the satisfaction functions of all the individ- 
uals are the same, and that within a group the amounts are divided 
equally. Therefore, an individual of the first group will possess amounts 
a1/N1, y1/Ni, and will produce «/N;. An individual of the second group 
will possess amounts a2/N2, y»/N2 and produce yo/Ns. Denoting now the 
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satisfaction functions of each individual by s; and s», correspondingly, we 
have 


x B 
$1= Aj log 37-+ Aa log 7 — oP; : (65) 
x 
se= Ay log =2+ As log 22 — > 20. (66) 
2 2 2, 


The total satisfactions of each group are given by 


x B 

Si Wisi= NiAy log 4+ Ni As log se — = to; (67) 
x J 

Sp = Nase = NoA1 log 7. + N2Ao log 7 — = 90. (68) 


To determine 21, *, y; and ye, we use the same procedure as before. For 
egoistic behavior we then obtain the following values for the optimum 
satisfactions after the exchange: 


Aja, A,;A2Q, No 


MEN A log ee ne 
es 08 edge dade a koe 8 (Aya) Na ok 2 Oe 
iy eign oe Ns to ea) (70) 
: "BCA, + 42) Na? 8 BCA Aa) 


For NV; = N2 = 1 these expressions reduce to expressions (39) and 
(40), as should be the case. 
The satisfactions per individual after exchange are: 


- Aja, AiA2d No A 
tim ules pay) 14? ea, pay et YD 
2 
A;A2Q, ee Aa (72) 


Stn Aa lOg oa As) N»2 Mera ta). 


As we see, both s; and s) depend on N;/N>2. What determines that ratio? 
For Ni/N2 = 0, we have 52m < Sim; for Ni/N2 = ©, Som > Sim. Hence 
there is a value of Ni/N>2 for which sim = Som. If Ni/N2 is less than that 

value, then Sim > Som, and the individuals of the second group will leave 
that group and join the first group, in order to increase their satisfaction. 
The reverse holds true if Ni/N2 is larger than the value for which 51, = 
~ Som. Hence WV. ,/N> will be maintained at that value, which makes sim = Sam. 
Putting 

Mi oa m Mi ied ui 9 3 

Tree he PAyae ’ ( ) 

No Ne 1 gy 


224 N. RASHEVSKY 


equating the right sides of (71) and (72), and rearranging, we find: 
A 
41 +1) log gi 41)= (4 41) log gt (74) 


Po Tg (Ay Ay) /(A FA) , (75) 


This gives 


From (75) we find 


Ai e (4,~ 42)/(A, FAQ) 


As 


y= oA = 
1 sere =A,)/(A,+A_) 


(76) 


As A;/A; > ©, — 1. The greater the satisfaction from the results of 
a given activity, the more individuals will choose that activity. 

It seems somewhat surprising that neither B nor a; and a, enter into the 
expression for 7. This is due to the particular choice of the satisfaction 
function. Exact solutions with other types of functions are very cumber- 
some. We may, however, give an example of a limiting case. 

Let the satisfaction due to the possession of an amount « of the first ob- 
ject of satisfaction be given by 


A, log(1+<), (77) 


while for the second object let the AEESESSETOD increase linearly, and be 


given by Agy. If 


B 
Aa (78) 


then an individual producing only the second object will put in all the 
available time T. If inequality (78) is reversed, he will not work at all. 
With these assumptions we find for the satisfaction functions s; and s2 


= a eee =f: oN ett EES 5 
Sy Ax log (1 =) As, ay Mi ) (79) 
Xo > EEG 
=AJA —- ———S- —— ———-—s 
52 = slog (1+ 52)+ 422 ds N3: (80) 


By multiplying s; and s, by N; and Ne, respectively, we find the satisfac- 
tions for each group as a whole. 


Now using the same procedure as before, which led to equations (34), 
(35), and (36), we find 


%1= 7% 5 %2= (1—7) xX; (81) 


A,\NoX%)(1—7) . 12. ANo%)(1—7) 
As(Ni tape! 2) laa Nee 


Ji 
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so that the total satisfactions for each group are: 


ce Xo A\NoX)(1—n) B : 
& Ny Ay log (1-452) + Web E03 (83) 
ne Xo _ ArNo%) (1 — 7») B 
So NeAy log (1+52)-+ day N give == Pe Vo- (84) 
We recall now that x» is determined by 
aS 
Ri 0, (85) 
whereas 
youl. (86) 
Introducing (83) into (85), we find: 
nN 4,4, +N, V72A2a2+ 4A ,a,B(1—7) (87) 


Not x= 


2B 


The other root, being negative, has no physical meaning. Since we are 
dealing here only with an illustration, we shall simplify expression (87) 
considerably by considering the case in which B is so small that the term 
in B under the radical can be dropped. This implies not too small an 7. 
We must have 


*Aja 
B.<i. (88) 
With this assumption, we find: 
NoAiain 
8 
BP ras (89) 


Now introducing (86) and (89) into (83) and (84), and dividing the lat- 
ter two by Ny and Ne, correspondingly, we find 


wa Ay log fet As 1A; (90) 
Aust AoT ss BT 91 
ee ae Nein (de Vs Acc 


Again n is determined by equating sim and S2m. This leads to the follow- 
ing equation 
Ar lon nT BT a $2) 
Go eo. Gada. 
Put 


See ala a eed a sare Rey oe (93) 
— =); Nowe: Aod, Be; B B’ ; 
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and again assume B so small that 
Saal Bez On (94) 


Equation (92) now gives 


2v + aB — Va?B?+ 4aby (95) 
ies 2v 


The other root is discarded as giving 7 > 1. 
Asy = A,/A,— ~, 1 — 1. This is seen by writing (95) as follows: 


mea (96 
ya 1S ae 
For very large values of v this becomes 
ap 
= 1—\— 97 
Tea a EAP 
Because of (93) this gives: 
T (Aza. —B) 
a nee ; 98 
Z : No Aid. ( ) 


Now n does contain A», as well as B,, No, and a2. The larger Wo, the total 
population, the more rapidly » approaches unity. 


ddd 


The results of the preceding section may be applied to the problem of 
the ratio of urban to rural population. Actually, both urban and rural 
populations each engage in a number of different activities. As a first 
rough approximation we may lump these activities in. two categories: 
agricultural and urban. The first produces mainly food. The second pro- 
duces other necessities, such as clothing, tools, etc. The ratio A1:/A» meas- 
ures the relative importance of the two. In early history the relative im- 
portance of food was most pronounced, Ai/A:2 was very large, and most 
of the population was agricultural. The ratio Ai/Ag, in general, also de- 
pends on climatic conditions. In a warm climate the need for such city 
products as clothing is much less than in a cold climate. Thus the ratio 
1— 7 of the urban to total population is also affected by climatic condi- 
tions at those early stages. 

As we have seen elsewhere (Rashevsky, 1952b) the value 1 — 7 deter- 
mines the rate of technological development. The increase in technology 
increases the value A». For example, in a modern society the need, even 
in the country, for such urban products as refrigerators, radios, etc. is very 
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great. Thus technological development, which increases with n’ = 1 — 7, 
favors in its turn the increase of 7’. 


Further study of these problems will be of importance for the dynamics 
of historical processes. 


The author is indebted to Dr. George Karreman for a discussion of the 
paper. 
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A mathematical theory of the process of the exchange of substances 
between the blood in the capillaries of a homogeneous tissue and the extra- 
cellular space, and between the extracellular space and the cells is devel- 
oped. An ideal geometry of the tissue is assumed, based to some extent on 
recent anatomical work concerning the functional distinction between 
two types of capillaries, the arteriolo-venular and the true capillaries. 
Equations are developed relating the concentration in the arterial blood 
to the mean capillary concentration, the concentration at the wall of the 
capillary in the extracellular space, and the average concentration in the 
extracellular space, and also relating the cellular concentration to the 
average extracellular concentration. The solutions of the equations are 
given for certain special cases and numerical results obtained. It is shown 
that the average extracellular concentration is a sensitive function of the 
permeability of the capillary wall and also is strongly influenced by the 
diffusion coefficient of the extracellular space. Furthermore, it is shown 
that the speed with which the average extracellular concentration ap- 
proaches the steady state is largely a function of the permeability of the 
capillary wall. 


I. Introduction. Most cells of organisms with closed circulatory systems 
must depend upon the system to bring to their local environments mate- 
rials necessary for their maintenance and to remove their diffusible prod- 
ucts from such environments. Since most cells are separated from the 
blood and lymph capillaries by a connective tissue matrix (Gersh, 1952), 
the efficiency of the supply and removal of substances (for any given ar- 
rangement of the capillaries) must depend partially on the characteristics 
of the ground substance of the connective tissue, and upon the amount of 
substance through which diffusion may occur, relative to the amount of 
fibrous and other material which is inert with respect to diffusion. 

The aim of the present work is to analyze the influence of the inter- 
stitial matrix and capillary permeability upon the exchange of substances 
between the cells and the blood. Many approximations and idealizations 
seem necessary to obtain results which may be compared with experiment. 
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However, as N. Rashevsky (1948) has pointed out in another connection, 
even if we knew exactly the geometry of a given tissue and knew precisely 
the blood flows, diffusion coefficients, etc. necessary for the calculation of 
the exchange relations, and could solve the exceedingly complex boundary 
value problem so presented, the results would still be of limited usefulness. 
This is true because the geometry, etc. of the solved problem would be 
special cases of an almost unlimited number of possible geometries. The 
analysis which is developed here assumes as known certain tissue parame- 
ters and certain data concerning blood flow in the capillaries, which are 
known only inexactly, or not at all. However, these data seem necessary 
for a detailed study of the exchange relations and are definitely data which 
could be obtained by experiment. The need for such data may stimulate 
work aimed at obtaining them. But, even in the absence of exact measure- 
ments, certain relations become evident and certain tentative conclusions 
may be drawn. 

In the following, a brief discussion is given of the labile character of the 
tissues and of the extracellular matrix, then blood-tissue fluid exchange is 
described. Next, consideration is given to some experimental and theo- 
retical work on blood-tissue molecular exchange. This is followed by a 
consideration of the anatomy of the capillary bed. On the basis of the pre- 
vious discussion, a tissue model is then presented and equations governing 
the blood-tissue exchanges are formulated. Some applications of the solu- 
tions are considered, and a discussion of various implications of the results 
is given. 

IT. The extracellular matrix and tissue changes. I. Gersh (1952) has re- 
cently discussed the characteristics of the ground substance of the connec- 
tive tissue. The ground substance is shown to be labile, and its plasticity 
is related in some cases to the action of certain connective tissue cells. 
Changes in the ground substance organization are shown to be related to 
tumor growth, regeneration, cyclic effects of endocrine glands, capillary 
permeability, immunological phenomena, and to calcification of cartilage, 
bone, and teeth. 

Other work (Ruzick, 1927; Hurst, 1933; Gellhorn and Regnier, 1936; 
Pearce, 1936; Lowry et al., 1942a, b, 1946; Horvath, 1946; Berg, 1947; 
Rhoades, 1948) indicates that the local environment of cells may change 
as a result of aging, tumor growth, regeneration and irradiation, and that 
these changes may in some cases be correlated with the rate at which cells 
may exchange substances with the blood. In this connection the work of 
L. Friedman and E. O. Kraemer (1930) on the diffusion of various sub- 
stances in gelatine gels may be mentioned. These workers found that the 
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magnitudes of the diffusion coefficients of several nonelectrolytes in a 
10% gelatine gel were only about one-half of their values in pure water. 

ITT. Blood-tissue fluid exchange. E. H. Starling (1895-96) formulated the 
hypothesis that the exchange of fluid between the blood and the inter- 
stitial spaces depends of the balance between the hydrostatic pressure 
excess of the blood, which tends to filter fluid into the tissues, and the col- 
loid osmotic pressure excess of the blood, which tends to filter fluid from 
the tissue into the blood. The mechanism depends on the relative imper- 
meability of the capillary wall to the blood colloids, and on the low hydro- 
static pressure of the extracellular fluid. Both points are confirmed ex- 
perimentally (Pappenheimer e¢ al., 1948, 1951; Field et al., 1932; McMas- 
ter, 1946). The work of E. M. Landis (1927, 1934, 1946), J. F. Danielli 
(1940), and J. R. Pappenheimer et a/. (Joc. cit.) has confirmed the details 
of Starling’s hypothesis. 

While Starling’s hypothesis, as discussed above, accounts for blood- 
tissue fluid exchange, and, together with a consideration of the lymphatic 
system, adequately accounts for the maintenance of blood volume, etc., 
the exchange of many substances between the blood and the extracellular 
spaces occurs so rapidly that only diffusion through the capillary wall 
seems to provide a mechanism which could give the required speed. 

IV. Blood-tissue molecular exchange. That the process of filtration 

through the capillary wall of substances in solution in the blood plasma is 
not rapid enough to account for the observed rate of transfer across the 
wall follows directly from the estimates of the maximum filtration rates 
observed. Pappenheimer et a/. (1951) calculate that the ratio of the rate of 
transfer by diffusion to the rate of transfer by filtration during very rapid 
filtration is about 116 for glucose, 83 for raffinose, and 27 for inulin, when 
the extracellular concentration is zero. The ratio becomes equal to unity 
only when the ratio of the concentration in the extracellular space to the 
concentration in the capillary is about 0.99 for glucose and raffinose and 
0.96 for inulin. Furthermore, the work of C. Hyman e al. (1952) provides 

- direct evidence that even when a tissue is becoming edematous the rate of 

capillary exchange is not measurably affected. 

Much experimental work has been devoted to the features of blood- 
tissue exchange (Dominguez e¢ al., 1935-36; Manery et al., 1939a, b, 
1941a, b; Hevesy and Jacobsen, 1940; Hahn and Hevesy, 1941; Gellhorn 
et al., 1944; Kruh¢gffer, 1946a, b; Flexner et a/., 1948; Morel and Marois, 
1949; Jones, 1951; Pappenheimer e¢ al., loc. cit.; Renkin, 1952). 

The experimental work mentioned above, with the exception of that of 
Pappenheimer ef a/. and of Renkin, was done with intact animals, and the 
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exchange relations observed are greatly complicated by the fact that vari- 
ous organs receive in any time interval very different amounts of blood per 
unit volume of the organ. Thus H. B. Jones (Joc. cit.) gives as the volume 
of blood per unit volume of tissue per minute received by the various or- 
gans of a man of approximately 70 kg. weight: thyroid, 3.0 to 10.0; kidney, 
5.0; brain, 0.50; marrow, 0.15; liver, 1.10; muscle, skin, fat, and bone, 
about 0.026. The work of Gellhorn et aJ. (1944) mentioned above also in- 
dicates that an injected substance is removed from the blood by the or- 
gans with the larger perfusion rates and restored to the blood as the sub- 
stance is gradually transferred to the organs with the lesser perfusion 
rates. It seems apparent from this data on the blood perfusion rates of 
various organs that the large differences in these rates may mask changes 
due to other differences between the tissues of various organs. In experi- 
ments in which the blood concentration of an injected substance is fol- 
lowed in an intact animal such seems to be the case. Work such as that 
of Pappenheimer ef al. (loc. cit.), which partially avoids these complica- 
tions by perfusing a relatively homogeneous tissue, may provide data re- 
flecting more clearly the influence of the capillary permeability and the 
extracellular matrix upon the blood-tissue exchange relations. 

The experimental data on blood-tissue exchange is usually described 
quite accurately by a sum of exponential terms, with two such terms be- 
ing sufficient in many cases, although occasionally one is adequate (Flex- 
ner ef al., loc. cit.), and sometimes as many as five terms are needed (Jones, 
loc. cit.). 

The time constants of the exponential terms have been interpreted in 
terms of total capillary wall area, and the over-all permeability coeffi- 
cients of the capillary wall, neglecting the distribution process in the tis- 
sue, and frequently neglecting the differences in the blood perfusion rates 
of the various tissues. Thus T. Teorell (1937a, b) reviewed the work of 
several earlier investigators and provided a theoretical model for blood- 
tissue exchange which does not distinguish between tissues with respect to 
their perfusion rates. 

However, the work of I. Bloch (1941, 1943) and of R. Smith and M. 
Morales (1944a, b) and of M. Morales and R. Smith (1944, 1945a, b, 1948) 
does take into account in a detailed fashion many important tissue param- 
eters and provides a rational interpretation of the experimentally deter- 
mined coefficients. The work of Smith and Morales particularly seems to 
provide a rationale for the investigation of the exchange of lipid soluble 


substances which takes into account the most pertinent physiological fac- 
tors. 
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The work of N. A. Michels (1936), B. W. Zweifach (1936-37, 1939, 
1940a, b, 1950) and R. Chambers and B. W. Zweifach (1940, 1944, 1946, 
1948, 1949) on the functional anatomy of the capillaries has made clear 
that the capillary bed of most tissues is composed of two types of capil- 
laries, whose functional significance must be taken into account in any 
complete theory of the blood-tissue exchange. 

V. The capillary bed. The work of Chambers and Zweifach (oc. cit.) has 
disclosed that, functionally, the capillary bed consists of central channels 
with side branches. Starting at the smallest arteriole, the central channel 
is composed of a metarteriole, with a discontinuous array of typical muscle 
cells along it, imparting vasomotion to the vessels which are about eight 
to fifteen microns in diameter. Next comes the proximal portion of the 
arteriolo-venular or a-v capillary with atypical muscle cells which cause 
changes in the diameter of the a-v capillary only under abnormal con- 
ditions. The distal portion of the a-v channel is invested with a coat of 
connective tissue. These a-v capillaries fuse to form nonmuscular venules 
about fifteen to twenty-five microns in diameter; these venules have a rela- 
tively heavy coat of investing connective tissue. Branching off the proxi- 
mal muscular portion of the central channels are precapillaries which have 
a proximal muscle-invested portion about twenty to thirty microns in di- 
ameter. The true capillaries anastomose freely and return to the distal por- 

tion of the a-v channels or to the nonmuscular venules. The angle at which 
_the first of the precapillaries branches from a given a-v capillary is more 
than ninety degrees. Later branches leave at successively smaller angles, 
and the true capillaries return to the a-v channels at a very slight angle. 
The junction of the a-v and precapillaries is marked by endothelial folds 
which act as sphincters, apparently controlling the flow into the true 
capillaries. The true capillaries show only passive changes in diameter, 
which result from changes in the blood flow through them. Figure 1, based 
on the work of Zweifach, gives a schematic representation of the capillary 
anatomy. 

The a-v channel normally has a continuous rapid blood flow through it, 
subject to recurrent vasomotion which consists of an irregular series of 
partial contractions of the metarteriole and the precapillary sphincters at 
intervals varying from thirty seconds to three minutes. Even when active 
flow through the true capillary is cut off, the capillary still has a fluid 
sin it. 

The angle of branching of the precapillaries from the central channel, 
and the relatively large diameter of the proximal portion of the precapil- 
lary, is apparently such as to maintain within the true capillary a more or 
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less uniformly low hydrostatic pressure; perhaps less, for the greater por- 
tion of its length, than the colloid osmotic pressure of the blood. 

Thus Zweifach and Chambers accept the general concept of Starling’s 
hypothesis—that filtration from the capillaries occurs where the hydro- 
static pressure exceeds the colloid osmotic pressure—but shift the locus of 
action so that filtration from the a-v capillary is believed to occur along 
most of its length, and absorption from the tissue spaces occurs along most 
of the length of the true capillary. 

The ratio of the number of true capillaries to the number of a-v capil- 
laries is a tissue characteristic and varies from less than 1 in the skin to 
more than 10 in the skeletal muscle (Zweifach, 1936-37). 
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Ficure 1. Schematic diagram of capillary bed 


The unselective permeability of the capillary wall in most regions of the 
body is adequately explained by the “pore” theory (Chambers and 
Zweifach, 1947; Zweifach, 1940b; Pappenheimer et al., loc. cit.). According 
to this, the endothelial cells of the capillary secrete a cement substance, 
which binds together the endothelial cells and constitutes a matrix be- 
tween them which is porous. The matrix permits the passage of low molec- 
ular weight substances with little hindrance. A decrease in the calcium ion 
content or the pH of the perfusing fluid softens the cement and enhances 
its dissipation, and thus results in an increase in the rate of filtration and 
diffusion from the capillary. In addition to the layer of endothelial cells, 
the capillary wall has an endocapillary layer, which is very thin and non- 
cellular, lining the inner surface. This may be a layer of adsorbed protein 
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(Zweifach, 1940b). A closely adhering, investing layer of argentophil 
fibrils, the pericapillary sheath, covers the outer surface of the capillary. 

P. Rous e¢ al. (1930, 1931), F. Smith e¢ al. (1931, 1932), and Zweifach 
(1940) have observed what they interpret as an increase in the capillary 
permeability from the arterial to the venous end. Danielli and Stock 
(1944) have disputed the interpretation but not the observation. How- 
ever, the reasoning of Danielli and Stock does not consider the functional 
anatomy of the capillary bed as outlined above, and hence must be ques- 
tioned. 

VI. A tissue model for blood-tissue exchange. On the basis of the foregoing 
discussion, an idealized model of a tissue will be constructed and the ex- 
change between the blood and the tissue in such a model will be consid- 
ered. The considerations are restricted to the exchange of substances which 
are lipid-insoluble. Substances which enter the red cells in appreciable 
amounts, which are appreciably adsorbed on blood proteins, or which un- 
dergo chemical changes at appreciable rates in either the blood or the extra- 
cellular spaces—all such substances require special considerations which 
would considerably complicate the equations to be developed and hence 
are excluded from the present treatment. Lipid-soluble substances may 
pass through the endothelial cells of the capillaries and also through the 
walls of the smaller arterioles and venules. Hence such substances are also 
excluded. The model applies only to a homogeneous tissue whose blood 
supply and ratio of extracellular to cellular phases are everywhere con- 
stant. The length of all capillaries, both a-v and true, will be taken as 
equal everywhere. The radii of the capillaries, and the number and loca- 
tion of the active capillaries will be assumed constant during the process 
considered. 

Since the capillaries form an irregular network with the arterial end of 
one capillary likely to be adjacent to the venous end of another, only vari- 
ations of concentrations in the extracellular spaces in a direction normal 
to the capillary axis are considered. 

_ Furthermore, the filtration rate along the capillary will be taken as con- 

stant, positive along the length of the a-v capillary and negative along the 
length of the true capillary. This is perhaps partially justified by the ob- 
servations of the gradient in capillary permeability mentioned above. 
Thus if the capillary becomes more permeable toward the venous end, the 
rate of filtration through the wall would be greater at the venous end, if 
the pressure were equal to that at the arterial end. But the pressure is 
lower, and hence, as an approximation, the filtration rate along the length 
of the capillary will be taken to be constant. 
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The number of parameters which must be taken into account is so 
large and the number of combinations of these occurring in the equations 
to be developed is so great that the Glossary of Symbols appended at the 
end of this work (see p. 258) may assist the reader in interpreting the 
equations. 

Distances along both the a-v and true capillaries will be designated by 
z. The length L is defined as the average functional length of the a-v 
capillaries, i.e., the length along which exchange between the blood and 
tissue may occur at an appreciable rate. At the arterial end s = 0 and at 
the venous end z =L . The true capillaries form an irregular anastomosing 
network in which the length of any one capillary is not readily defined di- 
rectly. But the ratio of the total length of all active true capillaries in a 
homogeneous tissue, i.e., those through which blood is flowing, to the total 
length of all a-v capillaries in the tissue could be determined. This ratio, 
which will be designated by A, and called the activity ratio, is assumed to 
be known. Then in our model we replace the actual distribution of capil- 
laries by arrays consisting of one a-v capillary, of length LZ, and A true 
capillaries also all of length Z, all in parallel. In the model, z = 0 where 
blood enters the true capillaries, and = L where blood flows out of them. 

We asssume that the mean velocity of the blood entering the a-v capil- 
laries in the homogeneous tissue is v;(0) and that entering the true capil- 
laries is 22(0). Since we consider only substances in solution in the blood 
plasma, it is convenient to express the blood concentrations in terms of 
concentrations in the red cell-free, protein-free fraction fi(z) of the total 
blood volume which is in the a-v capillaries at any point z along the length 
of the capillary. Let f2(z) be the corresponding fractions of the blood in the 
true capillary. Since the blood entering the two types of capillaries may 
have different numbers of red cells per unit volume, even f,(0) and f2.(0) 
may not be equal. Let 7; and 72 be the radii of the a-v and true capillaries, 
respectively. 

We will assume that the bulk of the material filtered into or out of the 
capillaries is water, so that the changes in the volume of an element of 
blood or perfusate are due solely to water filtration. The product of the 
filtration constant by the net differences between the hydrostatic and 
osmotic pressures acting across the wall of the a-v capillary will be de- 
noted by P;. The corresponding quantity for the true capillary will be 
denoted by Pz. It is assumed that P; > 0 and P, < 0 for 0 <2 
where L is the length of the capillary. This assumption is in accordance 
with the concept of Zweifach and Chambers, discussed earlier, that filtra- 
tion from the a-v capillary occurs along most of its length, and that ab- 
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sorption of fluid from the tissue spaces occurs along most of the length of 
the true capillaries. 


Then the material balance equation for any element of length dz and 
volume mrjdz along the a-v capillary is: 


mridfr(2)d2= [wrifi(s) (2) —arifi(et+dz) 4 (s+ds) (1) 
—2rr,Pidz) dt. 


The term on the left-hand side of the equation gives the accumulation of 
the fluid in the volume mrjdz, and the terms on the right are, in order, the 
amount of fluid entering the element, the amount leaving the venous end, 
and the amount filtering through the capillary wall. 

Simplifying (1) and using Taylor’s theorem, we have 


Ofte 6 Osi 1 2Py 
Che eT i © 


We will consider only the case where filtration is constant with respect to 
time. Then the solution of (2) is 


2P, 
fas) 0:(2) = fa(0) (0) —=s, (3) 


Thus, according to the assumptions made, an element of fluid moves more 
slowly as it approaches the venous end of the a-v capillary. Similarly, we 
find for the true capillary, 


ACA AOL: — Es (4) 


We now consider the variation along the length of the a-v capillary of 
a substance in the fluid of the blood. Let Ki(z,t) be its concentration in the 
cell-free, protein-free fluid. The amount entering an element of volume of 
an a-v capillary, rjdz, in time dt is 


mrifi() (2) Ki(s, t) dt. (5) 
The amount leaving the venous end of the element is 
arifi(atdz) n(s+dz) Ki(2+dz,t) dt. (6) 


The net amount leaving through the capillary wall is 


Qar,[(k, +4,P,) K, (2, ) — RE | dedi, (7) 
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where &, and hj are, respectively, the average permeability coefficients for 
exit and entrance of the substance through the capillary wall (these co- 
efficients may be different in the blood vessels of the brain and other spe- 
cial tissues), d; is the distribution ratio giving the ratio of the concentra- 
tion of the substance in the filtrate to the concentration in the protein- 
free plasma, and E; is the mean concentration of the substance in the 
extracellular fluid at the wall of the a-v capillary. The amount of the sub- 
stance accumulating in the element of volume is 


aridzd Ki (2, t). (8) 
Since we have asumed that no inactivation or chemical changes occur in 
the blood, the material balance requires 
arid 2d Ky (2, t) = [wrifi(s) 1 (2) Ki(, #) 
—rrifi(st+dz) y(s+dz) Ki(z+dz, t) | de (9) 
— [2ar, (k, +4,P,) K, (2, t) —2ar,RE|dedt. 


haha Hac 
Simplifying this, we have 


a2 2B = f(s) 01(2) Ki (2,4) — fi(z+dz) 1 (2+d2) 


5 : (10) 
X Ky (stds, 1) —7 (b+ d,P) K, (4, ) ds +7 Ei () ds. 


Upon using Taylor’s theorem and (3) we finally obtain 


0K, 0K, 
Ot 02 
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Essentially the same considerations apply to the true capillary, except 
that the filtration here is assumed to occur from the tissue into the capil- 
lary. Hence the net amount of the substance leaving through the wall of 
the true capillary in time df is 


2ar,[k,K, (2,1) — (ki —d,P,) Ei (t) |dazdt. (12) 


Here k, and ky are the average permeability coefficients for exit and en- 
trance of the material through the true capillary wall, d: is the distribution 
ratio of the filtrate, and Ej is the mean concentration outside of the wall 
of the true capillary. Thus, by the same reasoning that produced equation 
(11), we obtain an equation for K2(z,f). This is 


0K. 0 Ky 


; = Se (a) 8) 


2 2 
= Geol eet K, +> (hy — d,P,)E, (é) . (13) 
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We have neglected the radial concentration gradient in the capillary 
since it is likely to be very small. 

We now consider the concentration of the substance in the extracellular 
and cellular phases. We use the previously defined activity ratio A and 
assume that each a-v capillary is surrounded by a cylinder of tissue of 
average radius AR. The average extracellular concentration in this cylin- 
der will be called £, and the average intracellular concentration in the 
cells of this cylinder will be called C;. Similarly we assume that each of the 
true capillaries in our model is surrounded by a cylinder of tissue of aver- 
age radius R, with an average extracellular concentration £, and an aver- 
age intracellular concentration C2. If the number of a-v capillaries in an 
arbitrary volume V is NV, and if we know A, L, and N, then the radii AR 
and R are given by the requirement that 


V = [A°?R?+ AR?| w7NL+B=7RLNA(A+1) +B, (14) 


where B is the volume occupied by the blood vessels larger than capil- 
laries. The volume B is not likely to exceed about seven per cent of the 
total volume of any tissue. According to these definitions, AR + R is an 
approximate measure of the average distance between a-v and true capil- 
laries. 

We will also assume that the amount of material removed from the tis- 
sue by the lymphatic system is negligible; hence the considerations are 
restricted to tissues in which the lymphatic drainage during the time con- 
sidered is very small compared to the blood flow. 

The gradient of concentration of the material in the extracellular phase 
at the surface of the a-v capillary is approximately (cf. Rashevsky, 1948, 
pp. 15-23) 

2 (Ey —E1) (15) 

AR 
where AR/2 is the approximate distance between the a-v capillary and the 
site where the true concentration may be supposed to be about equal to 
the average concentration. Then at the wall of the a-v capillary we have 
E, —E;) 2D gl 

[(k,+4,P,) K,,,— RE ~ reer (16) 

where Km is the mean concentration in the protein-free plasma, D is the 
diffusion coefficient in the substance of the extracellular space, g is a factor 
given by the ratio of the cell-free area through which diffusion may occur 
to the total area (normal to the gradient of the concentration), and / is the 
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ratio of the shortest distance between two points in the tissue to the dis- 
tance which must be traveled between the two points without crossing any 
cell walls. The equation implies that the speed of diffusion in the extra- 
cellular space is very great relative to the speed of transport through the 
walls of the cells. Similarly, at the true capillary wall, we have 


2D gl 
[b,K,,— (ki —d,P,)E,) = (Bi—-E,) ©. (17) 
The amount of the substance flowing into the cells in the region of the 


a-v capillary in time dt is approximately 
a (A?R?— 72) [h,E, — hC\ 2d ade, (18) 


where /, and /y are permeability coefficients for entrance to and exit from 
the cells, s/v is the cell surface to volume ratio, and 0 is the fraction of the 
extravascular space occupied by cells. The expression implies that the 
speed with which the substance distributes itself inside of the cells is very 
great relative to the rate at which it crosses the cell wall. The amount of 
the substance flowing from the region surrounding the a-v capillary into 
the regions surrounding the true capillaries is given approximately by 

(B,—E:) 2D gl FEA dzdt. (19) 
Here 

2 (E, — Es) 
AR+R 


is the approximate gradient at the boundary of the region and 2 7A Rdz 
is the approximate area normal to the gradient through which flow must 
occur. 
The average net amount entering the region from the a-v capillary in 
time dt is 
2ur,[ (hk, +4,P,) K,,() —RE| ded. (20) 


The amount accumulating in the extracellular phase around the a-v 
capillary is 
TA he af d2dE, , (21) 


where a is the fraction of the extravascular volume not occupied by cells, 
and f, is the fraction of extracellular volume which is available for the dis- 


tribution of the substance, i.e., not occupied by fibers or other material 
through which diffusion cannot occur. 
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Then the material balance requirement, together with (18), (19), (20), 
and (21), gives 


A 4x DglAR 
w (AR? — 7?) af dzdE, = — (E, ech se dzdt 
t2mr | (ket dP.) UG (t) — kiE‘) d zdt (22) 
2P2 9 bs 
— 7 (A?R? — r?) (hE, —hC,} a dzdt. 
After cancellations and rearrangements, this becomes 
dE, _ 4D gl ( A 2171 (ki + d;P;) 
dt = af, (A?R?— r?) \A+ ;) @,-E,) "NERC ee: 
2riky bsh bsh ce 
VR p — OStci Sho 
af, (A2R? — py afe af _ 


Similarly, the amount of material flowing into the cells in the region of 
the true capillary in time dé is approximately 


a (R?— 72) (hE, — hl edt. (24) 


The amount of material flowing from the region of the true capillary to 
the region of the a-v capillary is approximately 


eae 25) 
(E, —E,;) 2D gl R R oad. ( 


Here 


is the approximate gradient at the boundary of the region and 27kdz is 
the approximate area normal to the gradient through which flow must 
occur. The average net amount entering the region from the true capillary 
in time dé is 

Der, (hk, Kp. (t) — (hl —d,P,) Ei) dedi. (26) 


The amount accumulating in the extracellular phase around the true 
capillary is 
a (R?— r2) afd 2d, . (27) 
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Then. since no material is inactivated in the region, (24), (25), (26), and 


(27) give 
7 L 
x (R— 92) af d2dk, = — GB) gal i ee 
eae 2 ARR 2 
(28) 
b 
— (k= d,P,) El) dzdt — x (R— 12) (hE, — h,C,] a adt. 
After cancellations and rearrangements, this becomes 
dE 4D gl 2 TR 
2 = (EOP Koki 
di ¢,.-) @_-®) AFD taf @- ma (1) 
(29) 
2r2(ko—dsP2)E,  bshy b she | 
Ge af (R73) af 2g a pameey 


Equations (23) and (29) imply that no net exchange of material occurs 
between one array of our model, consisting of an a-v capillary and its as- 
sociated true capillaries, and another such array. 

The concentration of the substance in the cells of the two regions must 
now be considered. The amounts of material flowing into the cells of the 
two regions are given by (18) and (24). The consumption of the material 
by the cells in the a-v region is assumed to be 


Git Qe (Co Goke (30) 


where g; is the consumption in—say, moles per unit cell volume per unit 
time—and q is a proportionality constant giving the change in the rate 
of consumption for a concentration different from some initial concentra- 
tion Co. Both gq; and g are negative for consumption. For production, 9» 
will be assumed to be zero and gq; positive. Then the amount of material 
produced or consumed in the a-vy region is given by 


a (AR? — 72) dzbf,lq, +9, (C,—C,) ldt, (31) 


where f; is the fraction of the cellular volume which is available for the 
distribution and production of the substance. The amount accumulating 
in the cells is given by 


a (A?R?— r2) dzbf,dC, . (32) 
Then (18), (31), and (32) give us 
m (AR? — 72) dzbf,dC, = 1 (A°R*— 7?) (hE, — foe dzdt 


(33) 
+7 (A2R? — rt) dzbf,lq,+ gq, (C,—-C,) ] dt. 
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Cancellation and rearrangement give us 


dC 5 
(Ce ni eile Oem (34) 


Similar considerations for the cells of the true capillary region give us 


dC ‘ 
EE 7 aps — hela + + 2 (C2) . (35) 
We may now calculate the concentration of the substance in the venous 
blood leaving the capillary array. We will assume that edema is not de- 
veloping. Hence the amount of fluid entering the array at any time must 
equal the amount leaving. Thus we have 


m7r77 (0) 0, (0) + Arif, (0) 7, (0) =ar2f, () v, Z) 


tAntii yy {L) 


From (36), using (3) and (4), we obtain a relation which will be used later: 


11Py 


A= : 
roP» 


(37) 
The amount of the substance leaving the a-v capillary in time dt is 


given by 
rr, (L) 0, (L) K,(L, t) dt, (38) 


and the amount leaving the true capillaries of the array is 


Anr2f, (L) 2, (L) K, (L, #) dt. (39) 


Then, after mixing the flows, the concentration of the protein-free venous 
plasma K, is given by 


ifr (L) 4 (L) Ki (L, t) +A rofe(L) 2 (L) Ka (L, 1) 


Nee ee 40 
AF, (0) 2, (0) +4727, (0) 0, (0) ey 


K, (t) = 


VII. The transformation of the equations. The problem of the blood-tissue 
exchange has now been formulated in equations (11), (13), (16), (17), 
(23), (29), (34), and (35). These equations, together with (3), (4), and 
(40), are sufficient to determine the concentration of a substance in the 
venous blood, at the walls of the a-v and true capillaries, and in the cells 
of the two regions, as well as the average extracellular concentration in the 
two regions—all in terms of the concentration in the arterial blood, the 
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blood velocity in the two types of capillaries, and the various tissue 


parameters. 
The solutions of the equations are most conveniently obtained by 


means of the Laplace transformation (Churchill, 1944). That is, we associ- 
ate with each function F(é) its Laplace transform f(p), defined by 


f(b) = fo em dt. 


We will let x/, x:, yi, w:, and wm; be the Laplace transforms of Ej, Ei, Ci, 
K;, and Kn;:, respectively, with 7 = 1, 2. 
The transform of (11) is 


dw,(z,p) , 2 
er += (P,— dP, — hy) w, (2, ?) 


(41) 


pu, — K,= —f, (2) 9, (2) 


re Ge 


where K» is the initial concentration in the capillary, assumed to be con- 
stant. Then 


dw 2 2k; 
f,(2) 9, (8) 2 = —[p-= ,— 4,P,—b,) |w, + f (6) + Ky. 4) 
Let us put 
pi = p= (Pi— dP, — ks) and f B 


Then, using (3), (42) becomes 


dw ia, az dz 
—Bw,+uei+K, f,(2)0,() f,(0) », (0) — wer (43) 
si 


The solution of (43) is 


w,(z, p) ae =f tae a fh poh 


Wy (P) B, — Ge) 


Fio % 

where wo(P) is the transform of the concentration in the protein-free plas- 
ma of the arterial blood, and fi0v10 = fi(0)2(0). The mean value of w(z,p) 
along the a-v capillary is given by 


W mi (p) Cr pete IG 1 2P12 \18:/2P, 
=7f O-eye as. 45) 


CRO Veet ati = T F109 %i¢ 
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Calculating the integral, we finally have, putting 


r 2 2 
M1 = gp lb —> (Pi diPs — ba) ] and eee 


? 
rif 10 V10 


Wl bt — MX = AG 1 
7 Sa Eee 1—(1- ity, 
wb, = MX, _ 1G i, Gi 4- 11) [ ( d,) ] o] 


or 
a _ Ko Wem bay aur fi, ple (1 —,)!™ 
eee. = M4) +e (1— eee 
(46) 
set) 


We will also need me 


w, (Ly p) =(w= 3°) (1a) of (1 aly +49. 47) 


The transform of (13) is 


d By 
pur, (2?) —Ky= — fy (2) (2) PP) 5% p _ py, (4, 9) 
(48 
2 / / 
+= (= d,P,) #4 (0) . 
Then 
d 2 
fe (8) 0, (2) = - [p-=,- k,) le, " 
2 
+7 = d,P,) %)(p) + Ky. 
Putting ; ; 
B,=P—,, P—f,).. Hp = (h,— 4.) 
and using (4), we have 
dws = dz ee 
8 ton ee K ie peeks ; (50) 
T2 


where foo%0 = f2(0)%(0). The solution of (50) is 


we (2, pei > ies 2Poz Vie cn 
w, (Pp) B— Me EG) 1 Fo 20 


Equation (51) implies that the concentration of the substance in ‘the 
protein-free plasma entering any true capillary is the same as that enter- 
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ing an a-v capillary. This is undoubtedly not correct, but may be approxi- 
mately so, since many true capillaries originate fairly close to the arterial 
end of the a-v capillaries. The mean value of w» along the true capillary is 


Wma (P) B2— wet, — Ko _ os B, (72/2 P2) 
w, (P) B, — u,, — K, rill Ce ee) oo 


Foo 29 


The right-hand side is equal to 


1 
Le ae 1+7, 
Xe ica sf (1— 2) 1%] , 
where 
AN bye gre a eee 
oy.2, ey | 2 - (P2 Ro) | and he = rates en ‘ 
Then 


a (eee et RON (oki lene ag 
5a (2) = (x, =) Ne (1 + 2) ) 


Ba (-— pibl= ta) Ko 
deo 


wy) yee 


We will also need 


w, (Ly ») =(w,- #2) a - e+ wi{1— (1—d)*] +32 


+3, 
The transforms of (16) and (17) are 


2Dgl 
[ (k,+4,P,) Oo en Rix] = (x, — x,) vhs 
and 
[tq — (8, d,P,) af] = (af — 2) OB 
Using the notations 
= l 
tes - ky — d2P2 AY 
y= Vo= aX 
Rit ae’ Ro : 
prea ey ie, 2Dgl 
AR(ki+4,P,)' ” Rie ? 


we have for (56) and (57) 
We et ey Oy 


ey eh = 
We —%a%_ p,X, = 0. 


(52) 


(53) 


(5 4) 


2 (55) 


(56) 


(57) 


(58) 
(59) 
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The transforms of (23), (29), (34), and (35), after the introduction of 
the Ey’s and C)’s as the initial values of the extracellular and cellular con- 
centrations, are: 


> 4D gl A ‘ 2r (ki + d;P)) 
i) a7, (AR 2) G 1) @ 7a! Gf, (A2R? —R2) 
; ; (60) 
27k , bshy She 
XO mi — Gf, (AR? — 7) See aie bay ne! 
Be 4D gl 2 roko 
Be 8 af (Oo) (Af 1) IY To, 
: ; (61) 
2 ro (Ro — doPo) , bshy Sho 
Ramer) 15 cio"? oaf,9 92! 
h — gx 
P11 C= <8 a. nae & + gyi; (62) 
h h goC 
py2—Co a es ; IO + goys. (63) 


The following notation will be used for the combinations of constants in 
(60) to (63): 


4g1D wae eer TD to\e 
ape oeen Cea! 50 Of, (R?— 1) Geek 


_ 2971 (ki + diP1) o ee EAI EE 
1G fg (A2R? — 7?) ’ 2 Gj agik? =f?) 


s Dats y 2 ro (ka — d2P2) | 
Mapa) aR, * 
Sisine athe See 
aaa ape Mies Pet (65) di— Wxko- 


Using these abbreviations in (60) to (63), we have 


, bf 
ph, — Ey Ee Eh Ps — Ba ERA e gee I (64) 


b 
px, —E, = — £,0, + £,0, +b 9 Hgts — 08+ oy 22 In, (65) 
py -— Co = 11 a eat S nee q2V1, (66) 
py. —Cy= 01 5 1 — oye + q2V2. (67) 
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From (66) and (67), we have 


afa 
CG 
ogee t ‘ot5 De (68) 
2 pto— Qe ! 


VA 


Introducing (68) into (64) and (65), we get, after some rearrangements, 


x, (p+é,+0,- 


r(P+e,+0 


If we let 


then (69) and (70) become 


Seer ts Oma vas 
(69) 
ae ( 
E 1 Co b fe (ete) =) 
7 | Um af, ) Ob oseD ; 
ae ear oa PW mo FHy%s 
(70) 
Chen «(o+2) 
“Eso o7, posto = vw 
0 ap por oy 
0102 
Oy = Poe epee yeepeare 
0192 
1 al 
8 ‘al page eek 8 alte) 
om ‘| Eo afa Presa 
Oy 1%, — 9B, +444) = B, (71) 
O,%, = £,%, — PW + WX, = B, - (72) 


Let us now introduce the symbols 


ee ae 1— (1— 2)! 
ees rad 


Ai (1 + y1) ‘ i= Xe (1 + 2) j 


2 = [ww — 3] 8+ Ko, as [o— 32] fabs + Ko. 
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Using this notation in (46) and (54), and collecting the equations (58), 
(59), (71), and (72) for convenience of reference, we have: 


2k) hws BZ (73) 
CANE G5 Waa) oe a Ae (74) 
yoo pe —w = 0, G55) 
Pte pve 0 eit), (76) 
VX +a,%,—&,x,—¢w,,=8,; C73) 
¥,%, — §,%,+a,*,—¢,w,,.=8,- (78) 


These equations, (73) to (78), together with (40), (47), (55), and (68), 
are necessary and sufficient to determine the venous concentration, Fi, Fp, 
Ci, C2, Hi, and £; in terms of the arterial concentration, the flow rates in 
the two types of capillaries, and the various tissue parameters. 

Using (75) and (76) in (73), (74), (77), and (78) to eliminate w,,; and 
Wm. we obtain | 


ee) oP Plate. pte — 2a, (79) 
Be ie le <2 )8— Py, 18, 0,85 = 2, (80) 
“1 (y,— oy.) +%, (a, +40) — &,%, = 8, (81) 
%5 (Wy — by%,) — Eo%, 1%, (a, +o, 05) = By - (82) 


VIII. The steady state. We may use equations (79) to (82) to get an es- 
timate of the relative importance of the various parameters in their in- 
fluence upon the steady state concentrations of a substance when the ar- 
terial concentration has been constant long enough for the steady state to 
be assumed. 

We use the theorem that if a limit exists for a function F(¢) as ¢ ap- 
proaches infinity, then pf(p) approaches the same limit as p approaches 
zero through positive values, where f(p) is the Laplace transform of F‘(#) 
and # is the transformation parameter (Doetsch, 1950, p. 458). Hence as- 
suming the existence of a steady state, which seems plausible on physical 
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grounds, we have for (79) to (82), where the various functions of 9, 1.e., 
a1, a2, Bo, Bi, Bo, 1, Y2, $1, and & must be evaluated at the limit as pesion 


E((m,(1— ¢,) 8] +2,8, 0, =Jim — 92, . (83) 
EL (1 — §) — 8,%] +£,B,0,=lim,— pZ, » (84) 
E\(,—9,) +E, (a, +4, 0,) —E,£, =Jim, 0B, » (85) 
E, (J —o»,) —E,t, +E, (a, +4,0,) =Jim,pB, - (86) 


As an illustration, we will consider the production of a substance in the 
cells at a constant rate g; moles per cubic centimeter per second. Hence we 
take g equal to zero. We will also take k; = ky = ki = ky, and n = fo; 
P,, d,P;, Po, and d2P: will be neglected relative to k;, since even during 
extremely rapid filtration into the tissue from the capillaries, ; is likely 
to be from twenty-five to two hundred times as large as P; or Py. The 
exact ratio depends on the molecular weight of the substance, the ratio 
being higher for the lower molecular weights (Pappenheimer ef a/., 1951). 
Since ); and ); are the ratios of the leakage through the capillary wall to 
the fluid flow into the capillary, we may put as a good approximation, 
when edema is not developing, 


1—)\y= ei and 1—)\,=e:. 


For while ; and d, may range from 0.001 to 0.2 (Danielli and Stock, Joc. 
cit.), the higher values occur only during the rapid development of edema, 
and we restrict ourselves to the case of no net filtration into or out of the 
blood. 

With these assumptions, we have the following values for the various 
coefficients in equations (83) to (86): 


ai = £15 a2 = £2; jim, pBo= 22 a; B= ty tay A 


ky k rifiod 
=— 1+ d +— ) — ne ) nae 22K /rsF 9%10) + 
V1 ] P; a6) 1 Py | a= MOG Rs (1— ) d 


1if29 V29 
Siergyy ee 


ea (2 Lky/r:F40%20)) 5 
lim — PZ pi (0 foe} Siro P10 (1 <a e~ 2 Lk,/rF 19% 10)) : 
po+o ; re ; 


jim pZo =e (0, co) 220740 (1 == e~ 2 Lk 3/7 F50%20)) é 
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Let us also assume that A = 1; this may be approximately the case in 
mammalian skeletal muscle under the experimental conditions used by 
Pappenheimer et al. (Joc. cit.), especially with the smaller perfusion rates. 
If we assume that fiov10 = foov20, then the coefficients of equations (83) and 
(84) become identical, as do those of (85) and (86), and Ei = Fi and 
£, = E,. This assumption is unlikely to be good, since according to the 
observations of Zweifach, fooveo < fiovio, but the lack of exact data on the 
relative flow rates, the fact that the assumption seems unlikely to affect 
the illustration qualitatively, and the simplicity gained seem to justify its 
use. Using these assumptions, equations (83) and (86) reduce, after some 
rearrangements, to 


7 Tif 10 V0 — p—(2 Lk, /74F 15% 19) eee 2Dgl 
Ey] Jee sae eee ie) Rh, 
(87) 
= — K (0;.% ee — e~ Lk ,/rF19%0)) 
1 
, R(R’ — r}) 
1 i a =- Diet cay ay ehee (88) 


Solving these for E; and Fy, we get 


B= K (0 ) ae bf =2=1| ~) - i](1 os e OB ihe) ) (89) 
Le ae 4 > f10 V0 ; 


T1 


E,=K (0, ~) + qbfs 
(90) 


x | =) — 1| 1 pitts (1 - eGbb yet) t 


We now need numerical estimates of the various parameters. Mam- 
malian striated muscle will be considered since this tissue is marked by 
large changes in activity. The average length of a capillary is about 0.7 
mm. and the average radius is about 5 x 10~* cm. (Rous ef al., 1930; 
Krogh, 1929). Jones (Joc. cit.) gives the average perfusion rate of muscle as 
0.01 to 0.030 volumes of blood per volume of tissue per minute. Pappen- 
heimer et ai. (1951) give an estimate of 70 cm.” per gm. of muscle as the 
capillary wall area. This estimate is much smaller than those given by 
Krogh (loc. cit.). Using equation (14) with B, the volume of blood vessels 
larger than capillaries taken as 7% of the total volume, we obtain 
R = 36.4 X 10-* cm. and 2N = 3.18 X 10° capillaries of both types per 
gm. Furthermore, we obtain v9 = 1.33 10-* cm. per sec. for the blood 
velocity corresponding to a perfusion rate of 0.02 per minute. Taking fio 
as 0.5, we have fiov10 = 0.67 X 10-* cm. per sec. Hence L/fror10 ~ 100. 


252 GEORGE W. SCHMIDT 


However, the estimate of 70 cm.? per gm. used in these calculations may 
be low. The value 140 cm. per gm. gives L/fiot10 ~ 50. This value will be 
used in all calculations and graphs. Naturally, v19 depends on the activity, 
but since more capillaries open up during activity, and since the values of 
both E, and E/ are more dependent upon the values assumed for R than 
those assumed for L/fiov10, we will keep this fixed and consider only varia- 
tions in R, Del, and fx. 

The values of Dg/ depend on the molecular weight of the substance con- 
sidered, on the state of the extracellular matrix, and on the geometry of 


Ol ' 10 100 
Ogl(cm.? secz')x107 


Ficure 2. Graph of Equation (90). Values of parameters are: L/fiovi0 = 50 sec.; 71 = 5 X 
10-4 cm.; k; = 1074 cm. per sec. 


the distribution of the cells in the tissue. The factor gi is likely to be about 
0.1, since g is about 0.15 and / is about 0.66. The values of &; are dependent 
on the molecular weight of the substance considered, and upon the state 
of the capillary wall, which may be altered as a result of various physio- 
logical changes. The data of Pappenheimer ef a/. (1951) permit an estimate 
of k; to be made for various molecules. Thus for glucose, k; ~ 10-4 cm. 
per sec., for sucrose, k; ~ 0.6 X 10-4, for raffinose k ~ 0.4 X 10-4, for 
inulin kj ~ 0.05 X10-*, and for NaCl, k ~ 1.6 X 1074. We will consider 
variations in Dgl from 10-° to 10-® cm2 per sec. and in & from 10-5 to 
10-* cm. per sec. 


Figure 2 shows the dependence of the ratio of E; — K(0, ~) to qf. 
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upon Dgl for four values of R. This ratio is much more dependent upon R 
than Dgl. Thus for a produced substance of about the molecular weight 
of glucose, the ratio is 400 for R = 15u, 1760 for R = 30u, 4000 for 
R = 45y, and 7300 for R = 60u. 

Figure 3 shows the ratio of E; — K(0, ~) to E{ — K(0, ~) asa func- 
tion of Dgl for various values of R. Rather large differences between the 
average concentration and the concentration at the capillary wall may 
exist for the smaller values of Dgl if K(0, ©) is small. Thus if K(0, ~) is 


2.4 


22 


2.0 


Ol ! 


0 
pgl (cm? secs!) x 107 


Ficure 3. Graph of the ratio of equation (90) to equation (91). Values of parameters are: 
L/frotio = 50 sec.; ry = 5 X 10-4 cm.; &; = 10 cm. per sec. 


zero, the ratio of E, to Ej for Dgl = 10-7 cm.’ per sec. is 1.038 for R = 15y, 
1.075 for R = 30u, 1.112 for R = 45y, and 1.150 for R = 60u. 

Turning now to a consideration of the average intracellular concentra- 
tion C,, we have from either equation (34) or (68) that in the steady state 
with constant production 4, 


vfs Ap (91) 


The variation in C, with constant q; is dependent on the average external 
concentration through the ratio of / to /. If we assume that /, and /y are 
of the same order of magnitude, and that the cell radius is about 107° cm., 
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then with f, ~ 0.5 and  ~ 107° cm. per sec. we have, roughly, that the 
concentration in the cell at any time is about seventeen times that concen- 
tration produced in one second in unit volume by a source of strength qu, 
plus the external concentration. The external concentration, as we have 
seen, is possibly several hundred to several thousand times the concentra- 
tion which would be produced in a unit volume in one second by a source 
of strength qu. 

Hence it appears likely that the intracellular concentration of a sub- 
stance produced in the cell may be influenced to a marked extent by the 
external concentration and thus by the diffusion coefficient of the extra- 
cellular material and the permeability coefficient of the capillary wall. 

IX. Transient conditions. In order to investigate the influence of the dif- 
fusion coefficient of the extracellular material, and the permeability coef- 
ficient of the capillary wall upon the speed with which a substance intro- 
duced into the blood distributes itself, we will again resort to some drastic 
approximations. We will assume that A = 1, ki = & = ki = kit, n = 1, 
fiotio = froteo, and that [i] < <1 and || < <1. Using the last as- 
sumptions, we find that 
1—(1—) 1 | 


i ee 
RES ane gCE a va 
and 
: cry phi CPs hates 

] — = 
ai Se tee Ma Gl +b 471) ca 

The set of equations (79) to (82) now reduces to 
—x'89, +28 p,= —Z,, (94) 
x (y,—$oy,) +, (a,—£,+¢,0,) =8,- (95) 


Let us consider the case in which the substance does not penetrate the 
cells, or penetrates them so slowly, relative to the speed of its distribution 
in the extracellular space, that we may assume that none of the material 
has entered the cells during the time required for a quasi-steady state to 
be attained in the extracellular space. Thus we assume that a substance is 
suddenly introduced into the blood in such a way that its arterial concen- 
tration remains constant. Then, if none of the substance was present ini- 
tially, we have 


Z,=—_; B,=0 (96) 


where K is the constant arterial concentration. 
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Solving equations (94) and (95) for a; and «1, using these assumptions, 
we have 


Whee K [ai — £1 + $191] , 
"1 P [ri (ar — £1 +191) + pr (Yi — 11) J’ sa 
Te ee K [yi — $1] (98) 


pl (a1 — £1 + ¢1p1) + pi (ti — $11) | 


By our assumptions we have o; ~ o2 ~ 0, and we also assume, as before, 
that |Pi| < < kiand|P;| < < &. Then the denominators of both equa- 
tions (97) and (98) reduce to 


Pp lip +rdip1— oii] - (99) 
We also have 
a1—§:+¢i1n=P+¢o1p1 (100) 
and 
¥1— $11 = — 1p. - (101) 
Using (99), (100), and (101) in (97) and (98), we get 
ess K (p + ¢1 pi) 
“1 blip +1d101 — orp? ey 
and 


K¢ip1 
Pp lnip +rid1 pr — $17?) © (103) 


x= 


Inverting these equations in p by routine methods of the Laplace 
transformation technique, we have, after some rearrangements, 


/— pee = to, D/a+ of s 
El Kj1 ppt er temaertonk (104) 


Ey = K{ = e— #(¢101)/(1+e,) } (105) 


A check on the first of these equations is provided by a comparison of 
‘(0) computed from equation (16) with £;(0) computed from (104). 
Both computations give 


jee (106) 
Le Dy 
We have from (105) 
dE; = Pipl e~ §(¢101)/A+e,) ( 107) 
dt Li 
and an 
Hy (0) _ iP 108 
Ene K ica (108) 
Figures 4 and 5 show 
ears 


1 p1 
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as a function of Dgl and fy, respectively, for four values of R. These graphs 
give the time at which £, attains approximately 63% of the capillary con- 
centration K, as well as the ratio of K to dE,(0)/dt. 

It may be noted that for a molecule of the weight of glucose, the time 
constant is 2.1 sec. for R = 15y, 10 sec. for R = 30y, 24 sec. for R = 45y, 
and 46 sec. for R = 60u. The time constant for the distribution process, 
for any fixed R, is markedly dependent on Dgl in the range 107° to 2 X 
10-7 cm2 per sec. and on & in the whole range considered. Thus for 


1000 


Rais 


et ' 10 100 
Ogl€m.2sec=') x 107 


Ficure 4. Graph of the reciprocal of the time constant of equation (105). Values of parame- 
ters are: afa = 0.1; 71 = 5 X 10-4 cm.; A: = 10-4 cm. per sec. 


molecules as large, or larger than, glucose the distribution time is very 
much dependent on the characteristics of the extracellular matrix. The 
permeability of the capillary wall greatly influences the distribution time 
for substances of any molecular weight. 

However, the approximation that \; and ); are negligible compared to 
unity renders the treatment of doubtful value for more than the very 
start of the process of distribution, for these approximations, of course 
imply that the concentration of the substance is eee the een 
along the length of the capillary and that the amount removed from the 
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blood is negligible compared to the amount flowing through the capillary 
in any given interval of time. 

X. Discussion and conclusion. The treatment of capillary exchange de- 
veloped above, while undoubtedly approximate and of limited validity, 
does indicate that both the diffusion coefficient of the extracellular matrix 
and the permeability coefficient of the capillary wall are important param- 
eters in the distribution of substances in tissue. Changes in the diffusion 
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Ficure 5. Graph of the reciprocal of the time constant of equation (105). Values of the 
parameters are: af, = 0.137 = 5 X 10-4 cm.; Dgl = 10 cm.? per sec. 


coefficient of the extracellular material are not so important as changes in 
the capillary permeability. Both parameters are subject to change, and 
changes in these parameters are followed by relatively large changes in 
both the steady state cellular and extracellular concentrations and changes 
in the initial rate at which a suddenly introduced substance distributes it- 
self. Since the metabolism of cells is undoubtedly influenced greatly by the 
concentrations within them of various substances, it appears that signifi- 
cant changes in cell metabolism may occur as a result of alterations in 
either of these parameters. . 
The consequence of a change in the ratio of the number of true capil- 
laries to the number of a-v capillaries has been investigated only in so far 
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as this may be correlated with changes in R. It is hoped that a more com- 
prehensive investigation of this problem will form the subject of a later 
paper. 

The writer wishes to acknowledge with gratitude his debt to the staff of 
the Committee on Mathematical Biology. Professor I. Opatowski has 
greatly helped the writer in developing the basic ideas which culminated 

in this work. Professor H. D. Landahl and Dr. George Karreman have, by 

helpful discussions, contributed considerably. The preparation of the 
manuscript has been the task of Mrs. Carlton Smith. The writer owes her 
a debt of gratitude which is difficult to express. 


GLOSSARY OF SYMBOLS 

A = the activity ratio, the ratio of the total length of all active true capillaries 
to the total length of all a-v capillaries in a homogeneous tissue, also the ratio 
of the number of active true capillaries to the number of a-v capillaries in the 
tissue model 

AR = radius of cylinder of tissue associated with each a-v capillary 

a = the ratio of the extracellular volume to the total extravascular volume 

b = the ratio of the cellular volume to the total extravascular volume 

Co = the initial intracellular concentration 

C, = the average intracellular concentration in the region of the a-v capillaries 
at any time t > 0 

Cy = the average intracellular concentration in the region of the true capillaries 
at any time ¢ > 0 

D = the diffusion coefficient in the substance of the extracellular space 

d, = distribution ratio, the ratio of the concentration in the filtrate from the 
a-v capillaries to the concentration in the protein-free plasma 

d, = similar to d;, but for the true capillaries 

Eo = initial extracellular concentration, assumed to be the same everywhere 

E, = average extracellular concentration in the region of the a-v capillaries at 
any time ¢ > 0 

Ey = average extracellular concentration in the region of the true capillaries at 
any time t > 0 

E; = extracellular concentration at the wall of the a-v capillary at time ¢ > 0 

E = extracellular concentration at the wall of the true capillary at time ¢t > 0 

fi = ratio of the volume of the cell-free, protein-free fluid of the blood to the 
total blood volume in the a-v capillary 


fio = fi(0) 
fy = similar to f, but for the true capillary 
Joo = f2(0) 


fa = ratio of the extracellular volume available for diffusion to the total extra- 
cellular volume 


fo = ratio of the cellular volume available for distribution of the substance to 
the total cell volume 
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g = ratio of the area normal to the concentration gradient in the extracellular 
space, through which diffusion may occur, to the total area 

h, = permeability coefficient of the cell wall for entrance of the substance into 
the cell 

hy = permeability coefficient of the cell wall for exit of the substance from the 
cell 

Ko = initial concentration of the substance in both the a-v and true capillaries, 
assumed to be constant 

A, = concentration of the substance in the a-v capillary, 0 < « < L, for the 
time ¢ > 0 

Kz = concentration of the substance in the true capillary, 0 < z< L, for 
time ¢t > 0 

Ki = mean concentration in the a-v capillary for any time ¢ > 0 

Kn2 = mean concentration in the true capillary for any time ¢ > 0 

K, = venous concentration of the cell-free, protein-free fluid of the blood at 
any time ¢ > 0 

k; = permeability coefficient of the a-v capillary for exit of substance 

ky = permeability coefficient of the a-v capillary for entrance of the substance 

k, = similar to ki, but for the true capillary 

k, = similar to kj, but for the true capillary 

LI = average length of the a-v and true capillaries, assumed to be equal 

= the ratio of the shortest distance between two points in the tissue to the 

average distance between the two points which must be traveled while not 
going through any cell 

N = total number of a-v capillaries in an arbitrary volume of a homogeneous 
tissue 

P, = product of the filtration coefficient by the net pressure acting across the 
a-v capillary wall, assumed to be constant along the length of the capillary 
and greater than zero 

P, = product of the filtration coefficient by the net pressure acting across the 
true capillary wall, assumed to be constant along the length of the capillary 
and less than zero 

p = Laplace transformation parameter with respect to ¢ 

gi = rate of production per unit volume, or rate of consumption per unit volume 
when concentration in cell is Co 

g2 = rate of change of consumption with respect to the cellular concentration 

radius of cylinder of tissue associated with each true capillary 

7, = radius of the a-v capillary 

7. = radius of the true capillary 

s = average surface area of any cell of a homogeneous tissue 


es) 
I 


¢ = time 

vy = average cell volume 

v, = velocity of the blood in the a-v capillary 
velocity of the blood in the true capillary 


V2 
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19 = 9,(0) 
v20 = v2(0) 


= Laplace transform of Ky 


hey SS 


Laplace transform of Ke 


Wm) = Laplace transform of Kit 
Wm2 = Laplace transform of Kime 


Wy = w(0,p) = w2(0,2) 
= Laplace transform of Fy 
«| = Laplace transform of E; 


22 = Laplace transform of EF» 
«, = Laplace transform of E; 
y, = Laplace transform of C; 


yo = Laplace transform of C2 
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z = distance along either a-v or true capillary, measured from arterial end 
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The distribution of membrane current in three models of nerve, when a 
longitudinal, linearly increasing current is applied, is derived. For the 
simple core conductor model it is shown that, if the region over which such 
a current is applied is large compared to the space constant of the model, 
the membrane current in the mid-portion of the region is a constant, inde- 
pendent of the distance, the time following the application of the current, 
and the impedance of the membrane. The effect of nonlinear membrane 
electrical properties is discussed. It is further shown that these conclusions 
apply equally to the case in which the simple model is surrounded by an- 
other concentric sheath (the double cable model). In this case the im- 
pedance of the sheath does not influence the membrane current in the 
mid-polar region. Finally, it is shown that the form of the solution for the 
saltatory model, for this type of applied current, is identical with that for 
the simple model. 


Most methods of causing electrical currents to flow radially across the 
membrane of nerve for purposes of studying the effects of such “‘polariza- 
tion” result in large spatial gradients of membrane current density in the 
vicinity of the electrodes. This difficulty has been overcome in some cases 
by the use of internal electrodes in nerve fibers which are large enough to 
accommodate them (Hodgkin and Huxley, 1939; Curtis and Cole, 1940; 
Marmont, 1949). A method of producing a truly wide polarizing electrode 
applicable to small nerve fibers was sought. This paper presents the math- 
ematical analyses associated with the successful outcome of that search. 

The necessary form of externally applied current such that the mem- 
brane current in the simple core conductor model is essentially constant 
over a large region has been known for a long time, viz., that the total 
longitudinal applied current should increase linearly with distance along 
the model. An approximation to such a distribution of applied current: 


* Present address: Department of Physiology, University of Chicago, Chicago, Illinois. 
265 


266 ROBERT E. TAYLOR 


was achieved by workers in M. Cremer’s laboratory as early as 1922 
(Cremer, 1929, p. 274) in an effort to determine an absolute constant of 
excitation. 
This type of applied current has been used by the author to investigate 
the effect of polarization on the conduction velocity of frog nerve. A pre- 
liminary report of this work has appeared (Taylor, 1950b) and the details 
of the experimental procedures and results will be reported elsewhere. The 


FicurE 1. Schematic diagram of filter paper electrode used experimentally. p-p: polarizing 
electrodes, s-s: stimulating electrodes, r-r: recording electrodes. Lower graph shows ideal cur- 


rent density in vicinity of nerve, which is closely approached by this electrode shape. [Com- 
pare equation (19)]. 


author has also had considerable success in using this method of applying 
current to produce local, nonpropagated electrically induced contractures 
in large regions of frog sartorius muscle (Taylor, 1951). 

Briefly, the electrode is a triangular shaped piece of filter paper (for 
nerve) thoroughly wet with Ringer’s solution and mounted in a moist 
chamber. Voltage is applied at the two acute corners as indicated in Figure 
1. The resulting current which flows along the lower edge of the filter 
paper, where the nerve is placed, increases in density in a linear manner 
from zero at the right-angled corner. 

As will become clear in the analysis to follow, the current which will 


MEMBRANE CURRENT IN NERVE 267 


flow through the membrane of a fiber in a nerve so placed will be from out- 
side inward if the direction of the current J and its first derivative dl /dx 
are the same, as in the illustration, and from inside outward otherwise. It 
will be assumed that the current density in the region of the nerve is pro- 
portional to the total longitudinal current flowing in some immediate 
vicinity of the fiber under consideration. It is shown in Part II that the 
presence of the resistive connective tissue sheath surrounding the nerve 
trunk does not invalidate this assumption. 

The model for the nerve which will be first considered is the usual core 
conductor model with an applied longitudinal current which increases 
linearly from zero at x = 0. It would be quite possible to solve the prob- 
lem for the exact type of applied current as indicated in Figure 1, but 


OUTER FLUID 


———— 


Ficure 2. Geometry of continuous core conductor model considered. Section of infinitely 
long cylinder. For nomenclature see text. 


since the electrical properties are assumed to be linear, it is more conveni- 
ent to utilize the principle of superposition, as will be seen later, so that it 
is only necessary to solve the problem for an applied current which in- 
creases from zero indefinitely. This will be done, following the derivation 
of a general integral solution which is applicable to any type of applied 


current. 
Part I. StupLeE Core Conpuctor MopEL 


A. Membrane current and voltage distribution in the simple core conductor 
model with linear, arbitrary membrane impedance in response to a constant, 
suddenly applied current of any form. Several treatments of the membrane 
voltage and current in the core conductor model, for various types of ap- 
plied currents, both for the transient and for the steady states, are to be 
- found in the literature (see Hodgkin and Rushton, 1946; Davis, in Lorente 
de N6, 1947 for references). The assumptions which will be made here are 
those introduced by Weber in 1873 (Weinberg, 1941), with regard to the 
geometry of the model. 
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Assumptions (see Fig. 2). 

1. The geometry of the system is an infinitely long cylinder with a con- 
ducting core, an external conducting layer, and a surface membrane with 
linear electrical characteristics. 

2. The potential is constant throughout the core and the external fluid 
at any point along the length of the model. 

3. The core and the external path are pure linear resistances. 

4, Any initial membrane voltage will remain constant and can thus be 
neglected during the analysis. 


Let: 


x be the distance along the length of the cylinder in centimeters 

t be the time in seconds 

p =s-+ jw be the complex variable associated with the unilateral Laplace 
transform; s > so, where so is large enough to insure the convergence of the 
integrals involved; 7 = ~/—1 

m = 2njf be the imaginary radian frequency associated with the bilateral 
Fourier transform 

V1 be the potential of the external fluid in volts 

V2 be the potential of the core in volts 

Vm = Vi— V2 be the potential across the surface membrane in volts 

Vn(«,p) be the Laplace transform of V,, with respect to time 

V(m,p) be the Fourier transform of V,,(x,p) with respect to distance, « 

r, be the resistance per centimeter of the external fluid in ohm/cm 

r2 be the resistance per centimeter of the core in ohm/cm 

%m be the resistance of one centimeter length of membrane in ohm cm 

v1 be the longitudinal current flowing in the external fluid in amperes 

72 be the longitudinal current flowing in the core in amperes 

I = 1, + 22 be the total longitudinal current in the model in amperes 

i, be the “polarizing current” per centimeter length flowing radially into the 
external fluid from without the model. Clearly 

eros 

ee 

tm be the current per centimeter length flowing inwardly across the membrane 
in amperes/cm 

Z(p) be the complex impedance of one centimeter length of membrane (see ex- 
plicit definition later) 

WHA ny: 


ils 


bp 


be the complex space constant 
; r 1/72 
X= lim A, = ( — 
poo ryt rs 
be the real space constant 
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A simple application of Ohm’s and Kirchoff’s laws to this model yields 
the following relations immediately: 


Differentiating V,, = V,; — V» with respect to x and substituting from 
the above, we have 
O Vin 
Ox 


= — fit + role , (1) 


and, since J = 7, + 2%, this can be written 


O) Wes 
Ox 


=12(71+ 12) ah Ty. (2) 


Differentiating a second time and substituting 


=e d 1 not 


gives the well-known differential equation for the core conductor model, 


viz., 
02 Vim 
Ox? 


ae (rit 12) tm = == 10s13 (3) 


This development has been written out to show clearly the origin of the 
continuity conditions. The potentials V; and V2 (and therefore V,,) must 
be continuous or an infinite current would flow through the nerve. The 
core current 7, must be continuous or an infinite current would flow across 
the membrane. Since 72 must be continuous, (2) states that the continuity 
of 0V,,/dx depends upon the continuity of J(). It is clear from physical 
considerations that if (a) is zero for large (positive or negative) x then 
V,, must approach zero for large (positive or negative) x. 

In order to solve equation (3) for any type of applied current I(x) 
(which, in general, may be a function of time), a relation between V,, and 
im in terms of the electrical characteristics of the membrane is needed. 
Equation (3) holds for any type of membrane (linear or nonlinear, func- 
tion of time and distance, etc.). If it is assumed, as it is here, that the mem- 
brane electrical properties are linear, i.e., are characterized by parameters 
which are not functions of time or current flow, and are not functions of 
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distance, the relation between V,, and z,, will always be a linear differen- 
tial or integral equation with respect to time only, with constant coeffi- 
cients. It may, however, be of infinite order (e.g.. polarization impedance). 

Under these conditions the complex impedance of the membrane will be 
defined in the following way: If we take the Laplace transform of the rela- 
tion between V,, and 7,, under discussion, we will obtain a linear algebraic 
equation in their transforms. On the assumption of all initial voltages and 
currents being zero, the ratio between the transform of V,, and that of tm 
is called the complex impedance Z(f). In terms of our nomenclature this 
can be written 


Lea (4) 


|x 


If in this expression the real part of the complex variable p is set equal to 
zero (p = s+ jw), the resulting expression Z(jw) is identical with the 
complex impedance commonly encountered in electrical engineering. The 
presence of the nonzero real part here is associated with the fact that dis- 
continuous driving functions are used (current applied suddenly) and 
amounts to the introduction of a convergence factor later on. 

With the introduction of the membrane electrical characteristics in this 
form, and the assumption that the external current 7, is suddenly applied, 
we can proceed to derive an integral solution for equation (3) in terms of 
an arbitrary Z(p) and z,(). The solution will then be carried one step fur- 
ther for the case that J = ax (where a is a constant) over a certain region, 
and some interesting properties of the solution for any membrane im- 
pedance will be pointed out. Finally, the complete solution for this type 
of applied current and a membrane with parallel resistance and capacity 
will be presented. 

The first step is to take the Laplace transform of equation (3) with re- 
spect to time. The Laplace transform of V,, with respect to time is defined 
as 


Vn, —) =Lil Vn(#,t)} =f” Vn (a, t) ePedl, (5) 
0 


and, similarly, the Laplace transform for in is im. 

By virtue of the known properties of the Laplace transform (see, e.g., 
Churchill, 1944), the transform of equation (3) for 7, suddenly applied (all 
initial currents and voltages zero) is 


d? Vim Py ITs 


mast. es (i+ lin = (6) 
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and substitution of 


Z(p) == 
gives 
d? le Fest [oc ( Tia ro) ety ra 7 
de® ~~ Z@) _ SS 
Since the complex space constant is defined as 
/ 
a= (Se 
Ag mae 
equation (7) becomes 
Ae es etary tm 
dx? — Wn = Ts . (8) 


The Fourier transform of V»(%,p) with respect to x is 
V(m, p) =Fe{ Vn (*, p) } =f, Vealeppyier de, 


and, similarly, the Fourier transform of 7, is 7. Taking the Fourier trans- 
form of equation (8) gives us 


WE 
Nm V— yaa, (9) 
Solving this for V: 


— ria 
es SEE 10 
Y= xm? — 1) De 
Equation (10) is the product of two transforms: 7,, the inverse of which is 
tl», and 


2 
= 11. 


the inverse of which is (Campbell and Foster, 1942, pair #444) 
Tire ene (12) 
2p 


The inverse of the product of two transforms can be expressed in terms of 
the inverses of the individual transforms by the application of the real con- 
volution theorem (see, e.g., Churchill, 1944). This theorem states that if 
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f(a), g(x), and h(x) have the respective transforms F(m),G(m), and H(m), 
and if F(m) = G(m)H(m), then 


AGE) = fo g(a) h(w—2)da= fg (x—2) h(a) de. 


Applying this theorem to equation (10), the inverse Vn(x,p) becomes 
immediately 


Vn (x, p) = oe “ip (2) exp (- 2a (13) 
Since by definition 
dae Vie 
PORES)” 


we can write down 


2 T1 
Lp eeP) ee a fie a) exp ( 


The complete solutions for the membrane current and voltage are then the 
inverse Laplace transforms of equations (13) and (14). Written out for- 
mally these are 


Vn (%,#) =L3{V,, (4, p)} = 


(14) 


ee 
ei Vw (5. 2) 82 mob OU Uh) 
T 


and 


8) +7@ 
in(%,) =L“in (2, P))=5- fim (@, p) dp, — t> 0. (16) 


Although equations (13), (14), (15), and (16) constitute the complete 
formal solution for the current through and the voltage across the mem- 
brane of the core conductor model as a function of time and space for any 
type of suddenly applied current and any (linear) membrane impedance, 
their explicit expansion may be simple, difficult, or impossible depending 
mostly upon the type of membrane impedance assumed. It shall be shown 
in the next section that from these expressions, certain general properties 
of the solutions may be elucidated without explicit expansion. 

If one is interested only in the steady state solution, it is not necessary 
to evaluate the integrals in equations (15) and (16). A well-known theorem 
(see Churchill, 1944) states that if the function f() has the Laplace trans- 
form F(p), and if pF(p) is analytic on the imaginary axis and the right 
half plane, then 


jim f(t) = lim pF (p). 
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Applying this theorem to equations (13) and (14) gives the general steady 
state solutions as 


11! m 


sia c2)) tr ey ie (4) e lA dz | iy) 


Une (x, ) = 


ial ie 
Cea ie ) ele-2lAd gz, (18) 
The complex space constant \, disappears because its limit as p > 0 is 
the real space constant X. 
These equations are readily integrated for any of the commonly en- 
countered types of applied current. 


B. Solution for applied current which increases linearly with distance; 
linear arbitrary membrane impedance. In the experimental work of the 
author (cf. Fig. 1), the applied current J had the following form: 


O05 e105 
CHD O ie las 
1%) =<" G24, iy KOS Wa (19) 
05 H1,< x, 
GF ¢< Oforall x. 


While the solution for this type of current is readily obtained, a great 
deal of space can be saved by the application of the principle of superpo- 
sition. This principle states that, if the electrical properties of any net- 
work are linear (characterized by parameters which are independent of 
time and current flow), and if the voltages V, and V, at any point in the 
network result from the application of currents J, and J, through some 
element, then, if J, = I, + J, is applied in the same element, the result- 
ing voltage at the corresponding point is V. = Va + V». 

Since the membrane voltage and current which results from the appli- 
cation of current of the form defined above for x > xo is well known 
(Hodgkin and Rushton, 1946; Davis, in Lorente de Né, 1947), it can be 
seen that the only type of applied current which must be considered to 
get the total solution is of the form 


0, i) <A 


Rea ae,— 2 t6s05 (20) 
preee few 
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The fact that J(«) becomes infinite does not cause any embarrassment 
since it appears only as the derivative. 


Substituting the derivative of (20) into equations (13) and (14) and in- 
tegrating gives the Laplace transform of the complete solution for this 
part of the problem, viz., 


T1Nc a 2/X, < 

7 fae. <0, 

Vn(%, p) = (21) 
ry. —a/X, Ss 
Bey ] ee 

and 
71a x/d 

ADA peace [ € ° ) x < 0 ) 
2p (ri + 12) 

in (#, p) = make (22) 
es roe —a/d, - > 0 
PACTery see eee OO 


These equations contain a result of great interest. Consider the region 
in which x is large enough to cause the quantity |e~**-| to be small. (It 
must not be forgotten that \. is complex.) If the real part of x/d. be de- 
noted by Re(«/d.), then 


| ene] = ga Rela) 


so that, e.g., if Re(w/X.) is greater than 3, |e~*/*«| will be less than 0.05. 
For large values of Re(x/X.), equation (22) becomes 


. 1, a 
eco ayia es 23 
bm (25:2) PORE) (23) 
The inverse Laplace transform of (23) is a step function in time, viz., 
Or i eli. 
in (™, t) = (24) 
MELE. 5 jas OY 5 
(rit re) 


This tells us that, if we are far enough away from the point at which the 
applied current begins to increase linearly, and not too close to the point 
where it ceases to increase, the membrane current is not a function of the 
membrane impedance, the time following the application of current, or 
the distance along the model. It is also to be noticed that the membrane 
current is directly proportional to a, the gradient of the applied current. 
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The restriction that the value of x be large compared to the reciprocal 
of Re(1/.) means that these conclusions will obtain over a greater or 
lesser region depending upon the form of the membrane impedance, Z(p) 
and the frequency range considered. It is clear that if the membrane im- 
pedance is of the same form as that of a network of resistors and capacitors 
(or of polarization capacities, see Cole, 1949) that the region will be larger 
for the transient state (high frequencies, p large) than for the steady state 
(p zero). 

While the membrane current is essentially constant and independent 
of the properties of the membrane if the region over which the longitudinal 
linearly increasing current is long enough, the current flowing longitudi- 
nally in the core is a function of the membrane properties and is an in- 
creasing function of x. Since 

j= tis 
mas 


the core current is given by 
12 = if tAOee 
The steady state core current is found by substituting equation (22) 


into this expression, multiplying by #, and letting » approach zero. The 
result is easily seen to be 


T1ar 


at? xz /d < 
2 ryaoray ‘ Os 
Ane 8 sai —2z/r > 
PEA tear a oo 


In the experiments of the author the current was applied over a distance 
of five centimeters of frog sciatic-peroneal nerve (Taylor, 1950a, b). It was 
found that over the region from x = 0.5 to « = 3.5 cm., the velocity of 
conduction of the nervous impulse was independent of distance for any 
given amount of applied current up to block. With lambda about 0.3 cm. 
for this nerve, the core current must have increased over this region by a 
factor of about 


CR 
Sg 
so that it is hardly possible that the current in the core per se had any ef- 


fect on the excitable properties of the membrane. 


C. Complete solution for linearly increasing applied current; membrane 
with parallel resistance and capacity. If the membrane is considered to be a 
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resistance and capacity in parallel, the relation between the voltage Vn 
and the current 7, at any point « is 
_ Vn COVn (26) 
Ts, ot ’ 
where C is the capacity per unit length of membrane in farad/cm. Under 
the assumption that the initial voltage and current is zero, the Laplace 
transform of equation (26) is 


Cis p) =1"4¢p vas 


and the membrane impedance is thus 


tee 
EADS ae eae 

where ¢,, is the membrane time constant frmC. 
In this case, then, the complex space constant is related to the real 


space constant by 


(27) 


2 
v= Roepe 
Substitution of equation (28) into the relations (21) and (22) gives the 
Laplace transforms, the inverses of which constitute the complete solutions 
for the membrane current and voltage for an applied current which is zero 
for x less than zero and increases linearly with « for x greater than zero 
(gradient = a). Making these substitutions, we have 


(28) 


r12a 


oN Vie 
: 2p 1 tiop) ae a ci 
Me 2 | (29) 
and : 
> a oVi+ m 
ENC er Pea LAG eS 0, 
im(", p) = ee pas 
ra £3 me e2V1+tmp/d | ; x = On 


l2p(n+7D) 


The inverse Laplace transforms for the various terms in equations (29) 
and (30) can be found in the tables of G. A. Campbell and R. M. Foster 
(1942), pairs #825.5, #448, #819, and #415. All of the pairs of transforms 
listed by Campbell and Foster as Fourier integral pairs can be used as 
Laplace transforms where the restriction is found that ¢ > 0 (their g). In 
pairs 7825.5 and #819 it is necessary to allow their parameter gamma to 
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approach zero. This is not permissible if these pairs are to be used as 
Fourier transform pairs, but can be shown to be legitimate if they are 
used as Laplace transform pairs, as in this case. 

With the use of these pairs the inverses of equations (29) and (30) can 
be written down by inspection. Letting 


x t 
AS Pa 
and introducing the error function: 


Y) z 
ase (2) =] loc) =} iiS= e dw 
tic (2) Tis) if 


Tv 


[tabulated values of error functions can be found in Tables of Probability 
Functions (1942)], the complete solutions are 


ria _x == Xe 
STi | « eric ( VF — 5" vid) 
3 ae eFeric ( — VE -5~m) 


—2e~Terfc =) , Fa oly 
i (x, 1) = (31) 


aS [4c —e-T) — eerie ( VF +59 
ns e~¥erfe (— VE 5 
+ 2 e~ Terfc Cal ; 220), 


and 


Frew: a [ - e-*erie (VF 5m) 
+ eXerfe(— VF-z~a)], RSA 
im (%, t) = (32) 
Sl eo [4+ ererte (VF+ 55) 
— eerie (— VF +549), pean We 


Equations (31) and (32) are the complete solutions for the voltage 
across and the current through the membrane with parallel resistance and 
capacity, for the simple core conductor model, when the applied longi- 
tudinal current has the form given by (20). It isa simple matter of addi- 
tion to get the complete form of the solution—if the applied current is of 
the form given by (19) as in Figure 1—by the use of equations (31) and 


Ficurr 3. Membrane voltage (in addition to resting) graphed according to equation (31). 
In A, time, and in B, distance, is the parameter. Simple core conductor model with membrane 
of parallel resistance and capacity when current of the form J(«) = 0 for « < 0 and ax for 
x > 0 is suddenly applied. Curves are normalized to unity at ¢ and x infinite. 
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Ficure 4. Membrane current graphed according to equation (32). Model and conditions as 


in Figure 3. x = 0+ and X = 0— indicate values of membrane current as X approaches zero 
through positive and negative values respectively. 
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(32), and the solutions for applied longitudinal current in which a spatial 
step occurs (Hodgkin and Rushton, 1946; Davis, in Lorente de N6, 1947). 

Equations (31) and (32) are plotted in Figures 3 and 4 for certain values 
of X and T. In every graph, the equations have been normalized so that 
the coordinates are dimensionless. Figure 5 shows the initial and steady 
state core current under the same conditions, i.e., applied current of the 
form given by equation (20). 

For the case actually used experimentally by the author (Taylor, 1950a, 
b) the normalized membrane current is shown in Figure 6. The conditions 


2.0 


1.5 
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Ficure 5. Core current according to equation (25) for steady state (T = ) and from inte- 
gration of equation (32) for T = 0. Model and conditions as in Figure 3. Normalized to be 


dimensionless. 


used to calculate this curve are those of Figure 1 and equation (19) where 
lambda is taken to be 0.3 cm., x) = 17A, and a, = 20X. Only the steady 
state is shown, but, as pointed out in the last section, the length of the 
plateau of membrane current is at its minimum in the steady state (cf. 


Fig. 4). 


D. Nonlinear membrane impedance and some general remarks. Any non- 
linear electrical properties of the nerve membrane would of course modify 
the above results. A nonlinear resistive component of the membrane im- 
pedance is undoubtedly present, as has been shown to be the case for non- 
myelinated nerve (Cole and Curtis, 1940; Hodgkin, 1947; Hodgkin, Hux- 
ley, and Katz, 1949) and for muscle (Katz, 1948). A detailed analysis of 
the nonlinear case would be extremely difficult, even if the membrane 
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properties of frog nerve were accurately known. It is reasonable to assume 
that the nonlinear voltage-current relation of the myelinated nodal mem- 
brane is similar to that of crab nerve or muscle, viz., resistance increasing 
slightly for applied anodal currents and decreasing considerably for ap- 
plied cathodal currents. The effect of such a voltage-current relation on 
the curves plotted in Figures 3 to 6 would be to distort the spatial transi- 
ent. For the case of applied cathodal currents this would be a favorable 
distortion, the reduced membrane resistance in the plateau region would 
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Ficure 6. Normalized membrane current for the continuous core conductor model in the 
steady state under conditions approximately as used experimentally [see Figure 1 and equation 
(19)]. Lambda is assumed to be 0.3 cm. so that a = 17 indicates that the triangular portion 
of the filter paper electrode was 5.1 cm. long. This graph was constructed by the superposition 
of curves from equation (32) and the known solution for membrane current in the vicinity of a 
unit step in applied longitudinal current (thin electrode at « = x). 


decrease the space constant \, sharpening the spatial transient and 
lengthening the plateau. The magnitude of the membrane current in the 
plateau region would not be altered. 

Conversely, the distortion for anodal applied currents would broaden 
the spatial transient and shorten the plateau. The effect of large anodal 
currents on the membrane resistance of frog nerve is not known, but for 
the crab nerves studied by A. L. Hodgkin (1947) the resistance was quite 
linear up to applied anodal currents of six times threshold. If this were 
true for frog nerve, the effect would be negligible for anodal currents, 
while being desirable for cathodal. 


There are some interesting relations of a general nature which help to 
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clarify one’s thinking in connection with the above analysis. Equation (3) 
is essentially one of two simultaneous partial differential equations in 
which the dependent variables are the membrane voltage V,, and current 
tm, and the independent variables are time ¢ and distance x. The relation 
between V,, and i,, in terms of the membrane electrical properties is the 
other simultaneous equation. It can be readily shown with the aid of the 
Fourier and Laplace transforms that if the second of these equations, ex- 
pressing the impedance of the membrane, is a linear integral or differential 
equation in x alone, with constant coefficients but of any order, then one 
can make the following statement: Given the solution of equation (3) for 
some driving function (the collection of terms not involving the dependent 
variables, in this case —77p), the solution for the case in which the driving 
function is the derivative with respect to x of the original one will be the 
derivative with respect to « of the original solution, at any time ¢. Now n 
is a constant and 7, for the case of the usual bipolar method of applying 
current is the first derivative of the 7, for the case here considered. Thus 
the solutions for these two cases are related as a function and its first de- 
rivative with respect to x at any time ¢. Since it is clear that the membrane 
current in the mid-portion of the region between two widely separated 
narrow “‘polarizing”’ electrodes approaches zero, it is equally clear that a 
plateau should occur in the mid-portion of the region considered above. 

Notice further that the membrane current is the first derivative of the 
core current. Thus the core current in the case of bipolar current applica- 
tion has the same form as the membrane current in the case of longitudinal 
linearly increasing applied current. Since it is obvious (or easily shown) 
that the core current in the mid-portion of the region between two thin 
electrodes will be almost independent of distance and of the membrane im- 
pedance, so must the membrane current be in this case. 


Part Il. THE DouBLE CABLE MODEL 


Since the results of these analyses were used in experimenting upon the 
sciatic-peroneal nerve of the frog, which is a large trunk surrounded by a 
dense and resistive connective tissue sheath, the effect of this sheath upon 
the current distribution has been considered. That this epineural sheath 
has a high resistance has been generally accepted for a long time (Bishop 
et al., 1926) and has been pretty thoroughly demonstrated by C. Rashbass 
and W. A. H. Rushton (1949). The sheath was removed and its im- 
pedance measured by the author (1950b), and any doubt which remained 
concerning this point was effectively removed by F. Crescitelli (1951). 


282 ROBERT E. TAYLOR 


The term “double cable model’ was introduced by Rashbass and Rush- 
ton who worked out several important relations. The assumptions in- 
volved are identical with those for the simple core conductor model treat- 
ed above, with the addition of another concentric membrane, hereinafter 
called the ‘‘sheath,’’ meant to represent the epineural sheath. A general 
integral will be derived which is applicable to any type of applied current, 
giving the steady state membrane and sheath current distributions. This 
integral will be explicitly solved for the case of a longitudinal linearly in- 
creasing applied current as above [Fig. 1 and equation (20)]. The geometry 
considered is shown in Figure 7. 


Ficure 7. Geometry of double cable model considered. Section of infinitely long cylinder. 
For nomenclature see text. 


The nomenclature is the same with the exception that: 

IT=i+ut+ ie 

7, = dI/dx is the current flowing radially into the outside fluid in amperes per 
cm 

Vo is the potential of the outer fluid in volts 

ro 1s the resistance per centimeter of outside fluid 

2 1s the longitudinal current in the outside fluid in amperes 

r, is the resistance of one centimeter length of sheath in ohm cm 

V, = Vi — Vois the potential across the sheath in volts 


Pens ae | ee 
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A. General integral steady state solution for any form of applied current. 
Using Ohm’s and Kirchhoff’s laws one writes down immediately 


dae otis (j= 0, 1, 2) 
di, 
m ax?’ 
: diy 
Us tm ax’ 
ent 
Up Las, 
Vi=t6%s 5 
Vie tv Peaks 


Differentiating V, = Vo — V; and V,, = V; — V2 and substituting, we 
have 


dV, _ r,dig } , 
Sf a rae — Polo Tih, 
1Vm _ mim 

ae = te — rl + Pete 


Differentiating again, solving for d7,/d«? and d%i,,/dx?, substituting in 
d,.” and X,,” and rearranging, we have 


Gig _ te imi to 
dx? 2 r Te 
8 8 $ (33) 
Che Ps Os — tet 4 tm 
dx? r ae 
m m™m 


Solving equations (33) by the use of the bilateral Fourier transform will 
allow the convenient use of the real convolution later. As before, the 
Fourier transform of z,, will be defined as 


Im (mM) = Bais (x) e~m*dx , 


where m is the imaginary radian frequency 27j/; 7./—1. Similarly, the 
Fourier transforms of 7,() and 7,(x) are 7.(m) and i,(m) respectively. 
Taking the Fourier transform of equations (33)}yields 


8 & : & (34) 
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and, solving these algebraically for 7, and mm, we have 


aay 
Rect G MS 


(35) 
Tori ES 
reser (a: 


The zeros of the denominator are the expressions defined above as 1/); 


2 
and 1/\;. The partial fraction expansions of equations (35) are readily 
obtained in terms of these expressions 


Introducing the following constants for convenience 


fe 
a eee 
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all of which can be shown to be real and greater than zero for any possible 


combination of physical values of the parameters involved, the partial 
fraction expansion of (35) is 


a. (on welled can C3 


m? — I m2? — Mia 
«) Gr) 
i, ( yeaa = 5 Lee e (37) 
Im\M) = rn ae os 1 ap? Bs ES Up - 
aoe V oe ) 


For any given form of the applied current 


dl 
hy (2) rk 
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equations (37) express the Fourier transforms of the membrane and sheath 
current distributions. Equations (37) however are each the product of two 
transforms and one can apply the real convolution theorem to this case as 
was done above in the analysis of the simple core conductor model. The 
inverse Fourier transform of z,(m) is by definition 7,(x), and the inverses of 
the multipliers of 7, are quickly found with the information that the in- 
verse of a transform of the form k/(m? — D?) is (Campbell and Foster, 
pair #444) 
ES <p) en Dla, 


provided only that D is real and greater than zero, which it is in every case 
under consideration. 


Thus the inverses of equations (37) can be written down immediately as 


Pe Oe ars fides e z-#I/¥1 — dg cze—l#- 41/2] 7,( 2) dz, 
1 ee (38) 
in (@) = BE Tyger = yen I] ig) ds 
These are the integral expressions for the sheath and membrane current 
distributions for the double cable model for any type of applied current 
i,(“), so long as the Fourier transform of the expression for such current 
distribution exists. Equations (38) are readily integrated for any of the 
commonly encountered types of applied current. 


B. Membrane current distribution in double cable model for longitudinal 
linearly increasing applied current. Integration of equations (38) for this 
case is very simple. Assuming an applied current of the form given by 
equation (20), i.e., J = Ofor« less than zero andI = ax for x greater than 
zero, 1»(“) is a unit step at the origin times a, the gradient of the applied 
current. The only part of the solution which will interest us here is the 
value of the membrane current for large values of x, since this will be the 
membrane current in the mid-polar region (cf. Fig. 1) if the total fein 
is more than a few lambda long. 

Integration of the second of the equations (38) under these conditions 
gives us 


i, (x) = Sai e~#/Ms) — )2(2 — e7 #2) J. (39) 
The plateau membrane current is lim 7, (x) which is 
x—- co 
206 ry 


i, (©) = (= 22), 
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or, substituting the value of c from (36), 


hee 

. Qroririr2 

in (@) =SOG 
137m 


However, from the definition of \; and Aa, 


x22 ft m 


2 Cro ry ry) 


and by substituting this and rearranging, the plateau membrane current 


is seen to be 


a 


tm (eo) = (40) 


ltn(S+2 
ay AY 

Thus the plateau membrane current for the double cable model with 
longitudinal linearly increasing applied current does not depend upon the 
resistance of the membrane or of the sheath and is linearly related to a, 
the gradient of the applied current. 

Integration of the first of the equations (38) is likewise not difficult, and 
it is readily shown that the sheath current in the plateau region is also 
independent of the membrane or sheath resistance. 

Integration of the equations (38) for the usual type of bipolar current 
application is straightforward. The resultant form of the membrane cur- 
rent distribution has been shown by the author to fit the data of Rashbass 
and Rushton (1949) on frog nerve for the spread of excitability in the 
vicinity of a single pole very well (Taylor, 1950b). These curves also fit the 
data published by Lorente de Né6 (1947) for the electrotonic spread in frog 
nerve and for the injury potential in the vicinity of a sharp cut. It follows 
therefore that these data of Lorente de Né furnish no argument for the 
contention that longitudinal currents in the nerve core influence the dis- 
tribution of membrane current in any except a purely passive, linear man- 
ner. No well authenticated data on frog nerve have been found by the 
author which do not fit the predictions of the double cable model theory. 
For an excellent discussion of some of these points the papers of Rashbass 
and Rushton should be consulted. 


Part III. Toe SALTATORY MoDEL 
A great deal of evidence has appeared to support the suggestion of R. S. 


Lillie (1925) that processes concerned with excitation and conduction in 
myelinated nerve occur at the nodes of Ranvier. At the present time no 
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phenomena which are clearly in conflict with this view have been present- 
ed. In any event the resistance of the myelin sheath is so large in compari- 
son to the resistance of the membrane at the nodes (Huxley and Stampfli, 
1949) that in a consideration of the effects of externally applied currents 
the assumption that the myelin is a perfect insulator is clearly justified. 

The purpose of this section is to inquire into what effect this type of 
structure will have upon the distribution of membrane current when a 
longitudinal linearly increasing current is applied (Fig. 1). 

The assumptions made here are identical with those for the simple core 
conductor model considered in the first part with the addition that the 
myelin is a perfect insulator, broken only at equally spaced interval by 


if”) intl) 


io(n+!)—— 


NODE n+! 


node width zero 


Ficure 8. Geometry of saltatory model considered. Section of infinitely long cylinder. For 
nomenclature see text. 


nodes of infinitesimal length and finite resistance (see Fig. 8). This model 
is identical with that considered by N. Rashevsky (1948) and Offner, 
Weinberg, and Young (1940) with reference to theories of conduction. 
Rashevsky’s expression [1948, p. 133, eq. (17)] for the electrotonic spread 
ahead of the active region contains a space constant identical with that 
which is here derived. 

The nomenclature is the same as that used in the first part for the simple 
core conductor model except that the unit of length shall be “node” in- 
stead of centimeters. If the distance between nodes is designated by d, 
then x = nd where n is the ordinal number of the mth node counting from 
the origin of x. Also, i,(#) and #,(m) are the currents flowing longitudinally 
between the node and the node m + 1 in the external fluid and core re- 
spectively, 7», is the total resistance of any node, and 7; and 7 are the total 
resistances of external fluid and core between any two adjacent nodes. 
These relations are indicated in Figure 8. 
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Applying Ohm’s and Kirchhoff’s laws to this model, we have immedi- 


ately 


Vi(n) — Vi(n+1) = ni (n), (41) 
Vo(n) — Vo(n+1) = rote (m) (42) 
tn (N) = (n) —an—1); (43) 
and by definition 
Vin (n) = Vi(n) — Von), (44) 
I(n) =i (m) +i2(n), (45) 
Via th) =t as (46) 


From (41), (42), and (44) we can write 
Vin (nm) — Vm(n-+1) = rit (m) — rete (m), 


and substitution of (45) gives us 
Vin (nh) = Bae (1) = 7h) = Cire) a) (47) 
Similarly: 
Vim (2-1) — Vin 4-2) oy (a Dr 1s) 1 eee 


Subtracting (48) from (47), substituting from (43) and (46) and introduc- 
ing the constant 


@=14 25% (49) 


Tm 


we find 


[I (n) —I(n+1)]. (50) 
This is the equation to be solved under the condition that a longitudinal 
linearly increasing current is applied, i.e., 


A/a <a es 
Con = Oe 


I(n+1) —I(n) —} (51) 


The solution to equation (50) under the condition (51) has the form 
(Norlund, 1924) 


oe nied Cy C1 = Cn EMM se n<Q 
m NG oer n>0, (52) 
where e”: and e” are the roots of the equation 
(lee 5 Ey ad bey 8 (53) 


and ¢y, c2, cz, and c4 are constants. 
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Solving equation (53) gives the roots as 


er=o+ Ve?—-1, 
op Se rg a 
so the solution (52) can be written 


Sesleae acl ne es Le 54 
m | 63.emn + cse EER ee ti me (0. (54) 

Before proceeding to evaluate the constants in equation (54) two theo- 
rems will be presented. 

Definition: A system S will be said to obey the principle of superposi- 
tion with respect to a variable s under the conditions that: If the driving 
forces fi(s), fo(s), and f3(s) produce the results 7;(s), r2(s), and r3(s) respec- 
tively, then for any f; and f2 such that fi + fo = fs, it is also true that 7, + 
Ly SPAR 

Theorem 1: In any system which obeys the principle of superposition 
with respect to the variable s, if f(s) produces r(s) and if f(s) is an odd 
function of s, r(s) is also an odd function of s. 

Proof: We can immediately write down the following: 

GS) produces <7(.s); 
=f (— 5s) produces: 7 (5). 
J (= 5). produces =r:(—-s)- 


Adding the first and third: 
f Gs) + f U— 8) produces -7(s) +7(—s). 
Subtracting the second from the first: 
f(s) +f/(—s) produces zero . 


Clearly then, r(s) + 7(—s) = 0. 

Theorem 2: In any system which obeys the principle of superposition 
with respect to the variable s, if f(s) produces r(s) and f(s) is independent 
of s, r(s) is also independent of s. 

Proof: The proof of this theorem follows directly when it is shown that 
if f(s) produces r(s), then df(s)/ds produces dr(s)/ds, provided the deriva- 
tives exist. The truth of the latter is immediately obvious if the derivatives 
are written out in terms of their fundamental definitions. 

We can now return to equation (54). Application of the above two theo- 
rems tells us that the solution i,,(7) has the property that 7,,(7) — in(0) 
is an odd function of ”. This is true because if the quantity a7:/27,, were 
added to the driving function of equation (51) it would then be an odd 


290 ROBERT E. TAYLOR 


function of m. But then the solution would be an odd function of » and 
would differ from the solution of (51) by not more than an additive con- 
stant. Since from physical considerations we have the boundary con- 
dition that, as 2 becomes large negatively, 7, must approach zero, clearly 
im must be bounded for large positive 7. Thus c, = ¢3 = 0 in equation 

54). After setting these equal to zero the bound for positive 7 is seen to 
be na/(r: + 7) and the above-mentioned additive constant must be just 
one-half of this and it must be equal to the value of 7,, (0). 

Therefore 
ra 


er 2 (rit 12) 
and from the second part of (54) 


. 2 2 ies Sie 
tm OTC ey Oe ae 
Thus 
piet ee es 
#2 (n+ 12) 
and (54) can be written 
71a 
——_——— e™, < 
f 2 (11+ 12) i ; aS 
im (n) = (55) 
pe ee» ea, SS 
a Ce ce ee 


If the quantity 1/7, is formally identified with lambda for the simple 
core conductor model, the form of the membrane current in the two cases 
is the same [cf. limit of (22) times p as p approaches zero], and all of the 
general conclusions apply to both. Thus the myelinated structure does not 
interfere with the plateau which is predicted for the mid-polar region when 
a linearly increasing longitudinal current is applied, with the direct rela- 
tionship between the membrane current density in this region and the 
gradient of the applied current, nor with the fact that the membrane cur- 
rent in this region is independent of the membrane resistance. 

The relation between the lambda for the simple model and 1/2, how- 
ever, is not simple and direct. It is easily shown that 


while 
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where the quantities in the latter case refer to one centimeter length of 
nerve and in the former to the distance between nodes. These relations 
tell us that the membrane at a node does not have the same effect on mem- 
brane current spread as the same resistance spread out over the distance 
d between nodes, and that measurements made on myelinated nerve and 
interpreted in terms of the simple model will lead to results which are 
consistent, but not to the correct values of longitudinal resistances or of 
the resistance at nodes (e.g., Lorente de N6, 1947, pp. 479 ff.). 


SUMMARY 


It has been shown that if a steady longitudinal current which increases 
linearly with distance along the nerve is applied to any of the passive 
models considered a plateau occurs in the curve of membrane current 
versus distance (or node number). The value of the membrane current in 
the plateau region is not dependent upon the properties of the membrane 
as long as they are linear, and is directly proportional to the longitudinal 
current gradient. 

For the simple core conductor model with membrane composed of re- 
sistance and capacity in parallel the complete solution for a suddenly ap- 
plied current of this kind has been presented. It is seen that in this case 
the region over which the membrane current is plateau-like is greatest im- 
mediately following the application of current and least in the steady 
state. From the general integral solution it was pointed out that this is the 
case as long as the membrane is such that its current voltage characteris- 
tics are like those in any equivalent linear network containing no induct- 
ance. It is clear that this conclusion holds for the double cable model and 
the saltatory model as well. Thus the steady state analysis gives the mini- 
mum distance over which a plateau in the membrane current would be 
expected to occur. 

It was pointed out that in the plateau region of membrane current, the 
core current is increasing linearly, with distance and that experimental re- 
sults strongly indicate that it is the transverse membrane current which 
modifies the properties of the nerve and not at all the longitudinal current 
in the core. 

Reference was made to experiments by the author that tend to bear 
out the applicability of these results to myelinated fibers of frog nerve 
with intact epineural sheath, viz., conduction velocity is altered, by ap- 
plied currents of this type, equally throughout the region of the expected 


plateau. 
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